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Abstract 

The search for materials displaying a large magnetoelectric effect has occupied researchers for many decades. The 
rewards could include not only advanced electronics technologies, but also fundamental insights concerning the 
dielectric and magnetic properties of condensed matter. In this article, we focus on the magnetoelectric effect in 
transition metal oxides and review the manner in which first-principles calculations have helped guide the search 
for (and increasingly, predicted) new materials and shed light on the microscopic mechanisms responsible for 
magnetoelectric phenomena. 



1. The challenge of discovering materials with large 
magnetoelectric responses - First-principles the- 
ory to the rescue 

1.1. Introduction 

The promise of new, oxide-based electronic de- 
vices has been fueling the interest in and search for 
new multifunctional multiferroics (see Figures [TJa) and 
b), a promising class of materials that possess more 
than one "ferro" property such as ferroelectricity and 
(anti/ferro/ferri)magnetism. One of the biggest chal- 
lenges in multiferroics research, which crosses many 
sub-disciplines of materials physics and chemistry, is 
finding or designing so-called cross-coupled multifer- 
roics that operate at room- temperature. These are ma- 
terials that are simultaneously magnetic and ferroelec- 
tric but in which the polarization (magnetization) can 
be manipulated in a useful way with an applied mag- 
netic (electric) field, i.e., those multiferroics that display 
a magnetoelectric effect, as shown in Figure [TJc). 

Structurally and chemically complex oxides have 
highly tunable ground states and are thus strong can- 
didate materials in which to realize new and enhanced 
magnetoelectric phenomena. The rewards could include 
not only new electronic device technologies, but new in- 
sights and challenges to our current understanding of the 
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condensed matter sciences. The family of known oxides 
is very large, but there are also an enormous number 
of possible tertiary and quaternary materials that have 
yet to be identified, much less characterized. Recent 
advances in thin-film growth techniques mean that we 
can now also create materials that do not nominally ex- 
ist in nature. Identifying promising candidate materials 
by surveying this vast space of design variables using a 
trial-and-error approach would be a hopeless task. 

In the past five years, theorists have made rapid 
progress in the application of first-principles methods 
to the design of new multifunctional multiferroic mate- 
rials. First-principles density functional methods have 
proved a powerful tool for studying the properties of 
materials at the level of atoms and electrons, without 
the need for empirical input. As such, the role of theory 
not only is that of analysis and interpretation of known 
experiments, but - increasingly - its focus is on predic- 
tion H H H HI El in [71. This has lead to more pow- 
erful ways of thinking about materials discovery where 
often theory is leading the way to the experimental dis- 
covery. In this article, we review the manner in which 
first-principles studies have increased our understanding 
of the magnetoelectric effect in transition metal oxides. 
We discuss the development of novel theoretical meth- 
ods to calculate the magnetoelectric effect and the new 
insights gained from applying first-principles methods, 
all with a vision towards the rational design of new ma- 
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terials. 

The organization of this review is as follows. In the 
remainder of the Introduction we discuss briefly the re- 
cent shift in our thinking about magnetoelectric design 
and a general strategy that has emerged. In Section [2] 
we focus the discussion on the linear magnetoelectric 
efl'ect, not only the methodologies to calculate various 
contributions to it, but also the strategies to design ma- 
terials with large linear magnetoelectric coupling are 
presented. In Section |3] we discuss spin-lattice and 
spin-phonon couplings and explain how they can be 
exploited to create new multiferroics and obtain diver- 
gent magnetoelectric responses. Microscopic coupling 
mechanisms and the efl'ects of particular exchange- 
correlation functionals in calculations are also briefly 
discussed. Section|4]involves first principles approaches 
to artificial heterostructures, which are designed to dis- 
play multiferroic properties. 

Although our review is comprehensive, it cannot be 
exhaustive and we therefore decided to omit the topic of 
magnetically-induced ferroelectricity, which would in- 
clude discussion on materials such as TbMnOs. Many 
theoretical models have been proposed concerning the 
origin of magnetically induced ferroelectricity in anti- 
ferromagnetic oxides and materials have been designed. 
We refer the reader to the original papers and reviews 
and references therein. E [9| |TOl HU [El HS] [H 

Another important topic that we do not cover is the 
use of first-principles-based eff'ective Hamiltonian ap- 
proaches |15|. First-principles calculations have been 
used for some time to obtain parameters for eff'ective 
Hamiltonians, which are consecutively solved using 
Monte Carlo methods (see ifTSl [T6J for well-known ex- 
amples). Although the use of such methods for the study 
of the magnotoelectric eff'ect is relatively recent, it has 
been very successful. Some important examples include 
the calculation of not only linear but also higher-order 
magnetoelectric couplings in BiFeOs 1 17 | and the pre- 
diction of the magnetotoroidic eff'ect, where a curled 
electric field can switch the direction of the magneti- 
zation |18|. Again, we refer the reader to the origi- 
nal references for more information on this interesting 

topic (nilllllll. 

7.2. Rethinking materials design 

Prior to the development of first-principles ap- 
proaches to the magnetoelectric eff'ect, symmetry ar- 
guments and a thermodynamic requirement were the 
primary guides for identifying materials displaying a 
large eff'ect. The thermodynamic requirement states that 
the components of the linear magnetoelectric tensor, a, 



should satisfy EOl 

where x and k are the magnetic and dielectric suscepti- 
bilities respectively. This lead to the general belief that 
a large dielectric response, that is a large k or equiv- 
alently a large dielectric permittivity 6, is required for 
a large linear ME eff'ect, a. Since ferroelectrics tend 
to be materials that display a large dielectric response 
(albeit in the paraelectric phase), a common misconcep- 
tion developed that magnetic-ferroelectrics, i.e., multi- 
ferroics, should be the targets of investigation for large 
linear magnetoelectric eff'ects. By definition, multifer- 
roics have both a spontaneous polarization and magne- 
tization and therefore lack both space and time-reversal 
symmetries. In many cases, a linear ME efl'ect would 
be allowed in such materials by symmetry. In the past, 
most design strategies have therefore focused on com- 
bining ferroic properties with the hope that a measur- 
able ME response could be observed. Although this ap- 
proach taught us much about new ferroic phenomena, it 
has enjoyed only modest success as means to identify 
materials for which there is a strong coupling between 
the various responses. There are a few difl'erent reasons 
for this. First, note that the thermodynamic relationship 
expressed in Eq.[T]is not a sufficient condition for strong 
coupling, that is, a large e does not necessarily lead to 
large a, it is merely an upper bound. Furthermore, it has 
recently been pointed out that known ME materials have 
magnetoelectric responses that are far from pushing this 
limit |[2Ti . In short, just because symmetry allows an 
eff'ect, it does not follow that the eff'ect will be strong. 

With this in mind, the focus recently has shifted away 
from combining ferroic properties and moved in the 
challenging direction of understanding the fundamen- 
tal mechanisms and key materials parameters that facil- 
itate a strong cross-coupled response between the vari- 
ous ferroic orderings. In order to create targeted mag- 
netoelectric phenomena, new approaches are being de- 
veloped that begin with devising a mechanism that con- 
trols the interplay between the diverse microscopic de- 
grees of freedom prevalent in these complex oxide fer- 
roics. Researchers are then utilizing principles of crys- 
tal chemistry, solid-state physics, and symmetry to de- 
velop sets of design criteria to aid in the identification 
of a real material with the physics of this mechanism 
built-in from the bottom-up. Finally, first-principles 
techniques, such as density-functional theory, are being 
used to screen for candidate materials and to incorporate 
materials-specific information into the conceptual mod- 
els that ultimately guide the search for new materials 
showing unprecedented magnetoelectric responses. 
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Figure 1: (a) An example of a multiferroic are ferroelectric- 
ferromagnets, which have a spontaneous polarization and magneti- 
zation, (b) Multiferroics are multifunctional in that they respond in 
a useful way to more than one external perturbation, e..g, an electric 
field couples naturally to the electrical polarization, P, and a mag- 
netic field couples naturally to the magnetization, (c) The generalized 
magnetoelectric eff'ect: The Challenge is understanding and subse- 
quently exploiting a microscopic mechanism that couples the electric 
and magnetic orderings at any order in the electric, E, and magnetic, 
H fields, (d) Cartoon representing the idea of phase competition. On 
the boundary between antiferromagnetic - paraelectric and ferromag- 
netic - ferroelectric phases, susceptibilities have divergent behavior 
and hence even small fields can induce colossal magnetoelectric ef- 
fect. 



1.3. A general strategy for ''colossal" responses 

Before moving on to the next section of this review 
it is worth discussing a general strategy that has been 
exploited in many different contexts to create multifunc- 
tional materials with large responses (and will be a com- 
mon theme in much of what we later discuss): the idea 
of phase competition. It has long been recognized that 
the key to enhancing sensitivity to an external perturba- 
tion is to position a material in the vicinity of multiple 
overlapping phase transitions (Bob Newhamm's Turn- 
bull lecture on Smart Materials | 22 1 and Tokura's work 
on CMR manganites |23 1 are classic papers). The mor- 
photropic phase boundary in Pb(Zr,TiO)3 - and the bi- 
critical region in (La,Sr)Mn03 are examples of regions 
where such overlapping transitions lead to highly re- 
versible phase transformations easily traversed by an 
external field and subsequently to an increase in the 
piezoelectric and magnetoresistive performance of the 
material. The main point we wish to highlight is that al- 
though the colossal responses displayed by Pb(Zr,TiO)3 
and (La,Sr)Mn03 are categorically diff'erent phenom- 
ena, they share a common "universal" feature - the com- 
petition among the respective microscopic degrees of 
freedom in each material leads to competing ordered 
states of matter. 

In regards to the magnetoelectric eff'ect, in 2006 
it was argued [24, _25 | that the interplay between 
spins and polar lattice degrees of freedom can be ex- 
ploited to produce a region of overlapping paraelectric- 
to-ferroelectric and antiferromagnetic-to-ferromagnetic 
phase transitions. The important thing about this strat- 
egy is that that the competition between these diff'erent 
ferroic orders in this phase competition region of the 
phase diagram, as depicted in Figure [ij^)' could lead 
to eff'ectively "colossal" magneto-electric and magneto- 
capacitive responses, without the need for the micro- 
scopic magnetoelectric coupling itself to be large. Ad- 
ditionally, it opens a route to achieve a large, non-trivial 
ME response from microscopic mechanisms beyond the 
linear ME eff'ect, i.e., at all orders of coupling, which 
removes many, if not all of the symmetry restrictions 
imposed by the linear eff'ect. 

The basic idea behind these design strategies to cre- 
ate strongly-coupled magnetic-ferroelectrics may seem 
counter-intuitive when considering the past obsession 
with multiferroics. Here the starting point can be in 
fact an antiferromagnetic-paraelectric (AFM-PE) mate- 
rial, which is the phase most magnetically ordered insu- 
lators have as their ground state. The challenge is two- 
fold: first, identify such AFM-PE materials that have 
a ferroelectric-ferromagnetic metastable phase (or in a 
previously inaccessible region of the phase diagram); 



and second, identify an external control parameter that 
can be exploited to tune between the various phases. 
Here first-principles approaches have a distinct advan- 
tage at mapping out this hugh phase space that is diffi- 
cult, if not impossible, to access experimentally - First 
principles to the rescue! 

2. The Linear Magnetoelectric Effect 

A discussion on magnetoelectricity in transition 
metal oxides inevitably starts with the linear magne- 
toelectric (ME) effect, a topic that has received much 
attention over the years ||26l[27l. This linear eff'ect in- 
volves the induction of an electric polarization that is 
proportional to an external magnetic field H, and like- 
wise, a magnetization that is proportional to an exter- 
nal electric field E. Pierre Curie was the first to predict 
the existence of such an effect over one hundred years 
ago; 1281 the first experimental observations were made 
almost 50 years ago |29|. Phenomenologically, it is re- 
lated to the lowest order (bilinear) coupling between E 
and H in the free energy |30|: 

r = ro- \KijEiEj - \xijtiiHj - atjEtHj^ (2) 

where Latin indices represent spatial directions, k 3.ndx 
are the dielectric and magnetic susceptibilities, and a is 
the linear magnetoelectric tensor, a is an axial tensor of 
second rank that is subject to very strict symmetry con- 
straints 1 31 J: only 58 of the 122 magnetic point groups 
allow a nonzero linear magnetoelectric effect. In addi- 
tion, except for 2 of these groups, multiple components 
of a are set to zero by symmetry (a detailed analysis of 
the form of a for all magnetic point groups can be found 
in 1 31 1). Given these rather stringent symmetry restric- 
tions, the fact that the linear magnetoelectric effect was 
predicted and subsequently experimentally confirmed to 
exist in a number of materials is a significant accom- 
plishment. 

Despite these early achievements, there is no known 
material that displays a large linear magnetoelectric ef- 
fect at room temperature. However, in recent years new 
first-principles methodologies have made it possible to 
calculate the linear magnetoelectric tensor in a relatively 
straightforward manner and in some cases with little or 
no modification of existing first-principles codes. 

Theory plays an important role in the search for new 
magnetoelectric materials by guiding experimentalists 
in their investigations and by helping to screen promis- 
ing candidate materials. In this section we discuss new 
first-principles methods developed to calculate the dif- 
ferent microscopic contributions to the linear ME effect. 
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Figure 2: Different contributions to the linear magnetoelectric effect 
at zero temperature. 



displayed in Figure |2] and the insights gained from the 
application of these methods. 



2.7. Advances in calculating the linear magnetoelectric 
tensor from first principles 

The linear magnetoelectric tensor, a, can be written 
as the sum of three terms: 



(3) 



where a^^^^ is the purely electronic part, is the ionic 
contribution, and a^^^^^^ is the strain-mediated magneto- 
electric response tensor, as shown in Figure [2] The elec- 
tronic contribution arises from the change in the elec- 
tronic polarization under an external magnetic field, or 
equivalently from the change in magnetization due to 
the change in electronic wave function under an exter- 
nal electric field. The ions are assumed to be fixed in 
their positions, as shown in Fig. |3jb). This clamped ion 
response can be measured at high frequencies where the 
ions cannot respond fast enough to changes in the ex- 
ternal field. Both a^^^ and a^^^^^"^ are lattice mediated, 
where the former is related to the change in magnetiza- 
tion/polarization with respect to the internal ionic posi- 
tions (Fig. ^c)) and a^^^^^^ arises from changes in the 
unit cell shape under external fields and so is related 
to partial derivatives of the magnetization/polarization 
with respect to the unit cell vectors (Fig. [3jd)). Each of 
the three contributions to a can further be decomposed 
into spin and orbital magnetization parts. We start by 
reviewing the studies that restricted themselves to the 
spin contribution as practical matter. 
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Figure 3: Sketch for different contributions to spin part of the linear 
magnetoelectric tensor: (a) The high symmetry system under no ex- 
ternal field, (b) The clamped-ion response to an external field, a^^^^, 
where the electron wavefunctions have changed but not the ionic posi- 
tions, (c) The change in ionic positions due an external field, accom- 
panied by a change in the electron wavefunctions due to this change. 
This contribution corresponds to Qf^°". (d) Lattice vectors can change 
in response to an external field, giving rise to of^^™". 



2.1.1. Spin-only magnetization at T=0 

The development of first-principles methodologies to 
quantify a (or its various contributions) was started in 
2008 by Iniguez |32|. In this work, a^^^^ was assumed 
to be negligible and only a^^^ was considered. It was 
shown that a^^^ has a very elegant form in terms of the 
Born mode eff'ective charges, Z^, defined as the deriva- 
tive of polarization, P, with respect to a phonon-mode 
amplitude, u, i.e.. 



dUn 



(4) 



where Pi is the component of the polarization, Un is 
the amplitude of the n^^ normal mode and Qq is the unit 
cell volume. In the past, there was serious debate |[33]| 
on how to define atomic charges in a solid. The Born 
mode eff'ective charge (and more generally the Born 
Eff'ective Charge |[34ll ) is well-defined, experimentally 
measurable 03, and has played a central role in the 
theory of dielectrics. For example, the ionic contribu- 
tion to the dielectric tensor has a very simple expression 
in terms of these mode eff'ective charges 1361 : 



^ne 



(5) 



where the summation runs over the infrared- active 
modes, e is the elementary charge and ojn is the fre- 
quency of the n^^ mode. 

What Iniguez did is to define a magnetic analogue of 
the mode eff'ective charge as the derivative of magneti- 
zation with respect to mode amplitude: 



(6) 



where Mi is the i^^ component of magnetization. Using 
this definition and considering Eq. |2] it is possible to 
obtain an elegant formula for a^^". 



-I ^"IK 



(7) 



which has an aesthetically pleasing (and intuitive) form 
similar to Eq. |5] (Here, Kn are the eigenvalues of force 
constant matrix.) Both and can be calculated by 
finite diff'erences using standard DFT codes, so the ex- 
pression for the ionic contribution to the linear magne- 
toelectric eff'ect can now be calculated in a straightfor- 
ward manner for systems of interest. Although Iniguez 
calculated only the spin contribution to Of, this method 
can also be used to calculate the orbital contribution, as 
has been done recently. 1 37 1 

This formalism was extended to the strain-mediated 
magnetoelectric response a^^^^^^ in Refs. |38 | and |39|. 
This contribution to a comes from a combination of 
piezoelectric and piezomagnetic eff'ects and can be ex- 
pressed in terms of the piezoelectric tensor e, the piezo- 
magnetic tensor h, and the elastic tensor C as 



eC-^h. 



(8) 



«=i 



The intuitive meaning of this expression is that, for ex- 
ample, an external magnetic field leads to a strain pro- 
portional to C"^/z, which in turn leads to a piezoelectric 
polarization proportional to eC'^h. So the resultant cou- 
pling between the electric and magnetic fields is propor- 
tional to the strengths of both piezoelectric and piezo- 
magnetic couplings. Note that all three tensors on the 
right hand side of Equation[8]can be calculated from first 
principles by finite diff'erences. Hence, standard DFT 
codes can be used to calculate a^^^^^^ as well. The or- 
bital contribution to a^^^^^^ can be evaluated by this for- 
malism too, but no such calculation has been presented 
yet to the best of our knowledge. 

Bousquet, Spaldin and Delaney calculated a^^^^ from 
first-principles for the first time and challenged the as- 
sumption that the purely electronic magnetoelectric re- 
sponse is small compared to the ionic contribution |40|. 
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Although it is necessary to use a full vector potential 
to simulate an applied magnetic field properly, the au- 
thors chose a computationally easier approach and con- 
sidered only a Zeeman field as an approximation. The 
sole eff'ect of the Zeeman field is to split the spin-up 
and spin-down states of electrons. This approach re- 
quires only minimal modification of existing publicly 
available DFT codes. In Ref. |40|, this is done for the 
VASP code 1 41 1 and it was shown that a^^^^ is compara- 
ble to a^^^ in Cr203. 

As discussed so far, the linear magnetoelectric eff'ect 
requires an interaction between the charge and spin dis- 
tributions of the material, in other words, between the 
orbital and spin degrees of freedom of the electrons. 
This interaction arises from spin-orbit coupling (SOC), 
which is a relativistic eff'ect 1421 |43l. In the reference 
frame of the electron, the nucleus orbits around it and 
hence creates a magnetic field that is proportional to the 
orbital angular momentum of the electron. This leads 
to an energy term that is proportional to L • 5^, where 
L and S are the orbital and spin angular momenta of 
the electron. This SOC term connects the orbital and 
spin degrees of freedom and breaks the spin rotation 
symmetry for the symmetries of the lattice. Not only 
is SOC responsible for phenomena like magnetocrys- 
talline anisotropy and the Dzyaloshinskii - Moriya inter- 
action (discussed in appendix), there would be no linear 
ME eff'ect at T = without spin-orbit coupling. 

2.7.2. Spin-only magnetization at finite temperature 

If the temperature is not zero, a diff'erent contribu- 
tion to the linear ME eff'ect emerges. Mostovoy and 
co-workers studied the eff'ect of spin fluctuations in 
Cr203 at finite temperatures using a combination of 
DFT calculations and Monte Carlo simulations {A^. 
The groundstate of Cr203 has no electric polarization 
because it is forbidden by symmetry (the point-group 
operation "space inversion followed by time reversal" 
leaves the system intact). At nonzero temperatures, 
however, it is possible that spin fluctuations break the 
relevant symmetries and lead to an electric polarization. 
An external magnetic field couples to the spin fluctua- 
tions and hence to the electric polarization. This leads 
to an extra contribution to a that exists only at nonzero 
temperature. An important difl'erence of this contribu- 
tion from the previously described zero-temperature ef- 
fects is that it does not rely on spin-orbit coupling but 
originates solely from the Fermi statistics and nonrela- 
tivistic Heisenberg exchange interactions. In the well- 
studied material Cr203, Mostovoy, et ah, found that the 
linear magnetoelectric efl'ect originating from the rela- 
tivistic efl'ects is about an order of magnitude smaller 



than the finite temperature contribution B4l . This sug- 
gests that the search for large a materials should be ex- 
tended to non-zero temperatures as well. 



2.1.3. Orbital magnetization 

All of the studies mentioned in this section so far con- 
sider only the spin magnetization of electrons. How- 
ever, as charged particles in motion, electrons in a 
solid can also have an orbital magnetization and a cor- 
responding magnetoelectric response B31 l46l l47l l48l . 
The orbital contribution to magnetization had not been 
considered in first-principles calculations until the re- 
cent developments in the modern theory of magnetiza- 
tion 1491 ESI EH E2 1- The clamped-ion orbital magneto- 
electric response tensor can be cast in a form that is the 
sum of a Kubo part and a pseudoscalar Chern-Simons 
part, which is isotropic. The Chern-Simons contribu- 
tion can be written as, 

^'='^'^^ (9) 

where e is the elementary charge, h is Planck's constant 
and dij is the Kronecker delta. ^ is a dimensionless pa- 
rameter (defined modulo 2n), which originally appeared 
in the context of axion electrodynamics 1 53 1. Axions 
are hypothetical elementary particles which were intro- 
duced to explain charge-parity violation |54|. If they 
were to exist, the Maxwell equations would be modified 
in a way that would lead to a bilinear coupling between 
electric and magnetic fields. While they don't exist as 
elementary particles in free space, axions have attracted 
a lot of attention in the field of topological insulators 
1551 . It is an interesting twist of fate that the idea of 
axions, originally conceived in the discipline of high- 
energy physics, is finding widespread application in the 
field of condensed matter physics. 

Even though a^^ has a compact form, there are 
many subtle details involved in the calculation of the or- 
bital magnetization (and therefore a^^ ) of an extended 
solid, which can be circumvented by the use of Wan- 
nier functions (561 E21 instead of Bloch states in calcu- 
lations. Coh and co-workers performed a detailed theo- 
retical analysis of a^^ using this Wannier function ap- 
proach [58 1 . 

While we were preparing this review a pre-print by 
Scaramucci and co-workers 1371 appeared on the ion- 
mediated orbital contribution to a. The first calculation 
of the Kubo terms for a real material will also be pre- 
sented soon by Malashevich, et al. ll59l 
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2.2. The power of first-principles methodologies: Elu- 
cidating new design strategies for achieving a 
large linear magneto electric effect 

Developments in the understanding and calculation 
of different contributions to the linear magnetoelectric 
effect reviewed in the previous section have lead to ef- 
forts to rationally design materials displaying a large 
linear magnetoelectric effect. 

Here we will see examples of the general theme dis- 
cussed in the Introduction to produce large responses. 
Specifically, strategies were developed to design materi- 
als with a strong linear ME effect for materials in which 
a is in fact relatively small and even zero by symmetry 
in the bulk ground state. The novel idea was to identify 
an external control parameter to tune the system to the 
proximity of a structural phase transition. 

2.2.1. Enhancing the linear ME effect 

Wojdel and Iniguez used the idea of 'structural soft- 
ness' in their first-principles study of magnetoelectric- 
ity in strained BiFeOs 1 39 1. They defined a material as 
structurally soft if its force-constants matrix {K) or elas- 
tic tensor (C) had a vanishingly small eigenvalue. Not- 
ing that the lattice-mediated contributions to a depend 
on K and C, they argued that structurally soft materials 
should have large values of a. This can be seen from 
Eq. ([t]), which indicates that a^^^ diverges as ap- 
proaches zero. Because 6^^" also diverges, positioning 
a system close to a ferroelectric phase transition should 
result in the desired enhancement. Selecting BiFeOs as 
their test material, Wojdel and Iniguez calculated and 
compared the stability, piezoelectric and magnetoelec- 
tric responses of two ferroelectric phases, namely the 
rhombohedral (R) and tetragonal (T) phases, as a func- 
tion of epitaxial strain. They found that for the R phase 
in particular, large compressive strains (-6%) induced 
a large increase in a, as shown in Figure |4] The au- 
thors argued that the enhancement is associated with 
a phonon mode with a vanishing force-constant. Ide- 
ally, the soft mode associated with the enhancement of 
a should have a large mode effective electric charge 
and magnetic 'charge' Z^. As noted by Wojdel and 
Iniguez, this is not the case with BiFeOs, but a signifi- 
cant enhancement is observed nonetheless. 

2.2.2. Inducing the linear ME effect or weak ferromag- 
netism: Two sides of the same coin 

First-principles calculations also recently predicted 
a large magnetoelectric effect in the antiferromagnetic 
insulator CaMnOs |60', ^6T\. CaMnOs adopts the or- 
thorhombic Pnma structure in its ground state but first- 
principles calculations 1601 revealed the existence of a 




Epitaxal Strain 8 (%) 

Figure 4: First principles calculations of structural softness leading 
to an enhanced a) piezoelectric and b) magnetoelectric response in 
BiFeOs. Near the structural phase transition at large compressive bi- 
axial strain, the frozen-cell (light) and full lattice-mediated (dark) re- 
sponse are enhanced in both the R and T phases. From l39l . Copyright 
(2010) by the American Physical Society. 

(weak) zone-center ferroelectric instability in the cu- 
bic Pm3m phase, which is suppressed by the octahe- 
dral rotations in the ground state. Bhattacharjee and co- 
workers have shown however that a proper ferroelectric 
instability could be induced by tensile strain. Addition- 
ally, they showed that the ferroelectric instability is Mn- 
dominated, contrary to earlier arguments f62| that dis- 
placement of the B-site cation from its ideal Pm3m po- 
sition would not be favored in materials with partially 
filled d orbitals. Hence, Bhatacharjee, et al, demon- 
strated that magnetic ions could contribute to ferroelec- 
tric distortions, an important development in the search 
for multiferroics. 

In 2011, Bousquet and Spaldin [61 1 showed that 
strain can induce a linear magnetoelectric effect in 
CaMnOs by inducing ferroelectricity. In Pnma 
CaMnOs symmetry allows a small canting of the anti- 
ferromagnetic spins, that is, weak-ferromagnetism, but 
forbids a linear magnetoelectric effect (a modest effort 
can be used to show that in a material with a center of 
inversion, the linear magnetoelectric effect and weak- 
ferromagnetism are mutually exclusive; see the Ap- 
pendix). Using first-principles calculations, Bousquet 
and Spaldin showed that Pnma CaMnOs can be driven 
ferroelectric with strain. Since a ferroelectric distortion 
removes the symmetry that forbids the existence of the 
linear magnetoelectric effect in the paraelectric struc- 
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ture, the linear magnetoelectric tensor is proportional to 
the ferroelectric polarization and therefore changes in a 
can be quite large (note that since the ferroelectric dis- 
tortion may in fact change the type of collinear, antifer- 
romagnetic order, this conclusion of ferroelectrically- 
induced linear magnetoelectricity is only true if in the 
ferroelectric Pmcli structure, the collinear spin struc- 
ture remains that which allows a linear magnetoelec- 
tric effect; Bousquet and Spaldin explicitly confirmed 
this for CaMnOs). The symmetry analysis presented in 
Ref. 1 61 1 is valid for all Pnma ABX3 perovskites. Con- 
sidering that Pnma is the most common space group 
among the ABO3 perovskites, Bousquet and Spaldin 's 
approach represents a promising route to search for new 
magneotelectrics. 

In 2008, Fennie |I63] |64l formulated a set of struc- 
tural and magnetic criteria that must be satisfied in a 
particular material in order for a ferroelectric distortion 
to induce weak ferromagnetism, a idea first suggested 
back in 1977 by Fox and Scott |65|. Quite simply, 
the criteria are the complement of the ideas expressed 
by Bousquet and Spaldin in 2011. Fennie argued that 
the paraelectric-antiferromagnetic reference structure 
should allow a linear magnetoelectric eff'ect, which sub- 
sequently forbids weak ferromagnetism. Then, the fer- 
roelectric distortion could induce weak ferromagnetism 
and as such the induced magnetization would be propor- 
tional to the ferroelectric polarization. Using FeTiOs as 
the test case, Fennie then used first-principles calcula- 
tions to show that both criteria are satisfied in a series 
of R?>c ATiOs materials (A=Fe, Ni, Mn). Experiments 
subsequently verified that R?>c FeTiOs is both ferroelec- 
tric and a weak ferromagnet 166 1. These predictions 
sparked a search for other materials satisfying the same 
criteria, a number of which have also been recently syn- 
thesized [67 , 68 ,69J. 

The work of Fennie and that of Bousquet and Spaldin 
showed the utility of using a somewhat "boring" para- 
electric reference structure to search for "exciting" phe- 
nomena in ferroelectric distorted ground states. We 
point out that the symmetry in the ground state of 
the materials considered by the diff'erent authors al- 
lows ferroelectricity, weak-ferromagnetism and the lin- 
ear magnetoelectric eff'ect. However, qualitatively dif- 
ferent responses to applied strain occur in each material, 
which becomes clear when using a higher- symmetry 
reference- structure. For example, increasing the po- 
larization in strained CaMnOs would, to lowest order, 
increase the linear magnetoelectric eff'ect while not in- 
creasing the net magnetization, whereas the exact oppo- 
site would occur for FeTiOs (in the Appendix we pro- 
vide a simple Landau model to show this). 



We have discussed an example where a ferroelec- 
tric distortion induces the linear magnetoelectric eff'ect 
and an example where a ferroelectric distortion induces 
weak ferromagnetism. Can a single material display 
both eff'ects? What we mean is, can a ferroelectric dis- 
tortion induce both the linear magnetoelectric eff'ect and 
weak-ferromagnetism in a structure that displays nei- 
ther? In order for this to happen, the ferroelectric dis- 
tortion must change the translational symmetry of the 
material, which is not possible in a proper ferroelectric 
transition like that in strained CaMnOs and FeTiOs . It 
could occur however in an improper ferroelectric transi- 
tion. Recently, Benedek and Fennie fTOl have discov- 
ered such a material: the ^ = 2 Ruddlesden-Popper 
material Ca3Mn207. They showed that the ferroelec- 
tric distortion was driven by two diff'erent octahedral ro- 
tation patterns with diff'erent symmetries |70|. Specif- 
ically, there is a trilinear term in the free energy of the 
form, 

r = yPRiRi, (10) 

where P is the polarization and Ri and R2 are the 
octahedral rotations. The combined rotation pattern 
Ri © R2 is treated as a 'hybrid' distortion and is what 
drives the system into the polar state, that is, the 
polarization becomes non-zero only when both rota- 
tions condense fTTIl . For these reasons, Benedek and 
Fennie coined the term hybrid improper ferroelectric- 
ity to describe this mechanism, which was first dis- 
covered by Bousquet, Dawber and co-workers 1721 in 
SrTiOs/PbTiOs superlattices. In the case of Ca3Mn207, 
the octahedral rotations induce not only ferroelectric- 
ity, but also weak ferromagnetism and a linear magne- 
toelectric eff'ect. Rondinelli and Fennie subsequently 
elucidated the design criteria for the creation of new 
hybrid improper ferrroelec tries using non-polar Pnma 
perovskites as 'building blocks' 1731 . As mentioned 
previously, Pnma is the most commonly adopted space 
group among ABO3 perovskites 11741 , hence the study 
and synthesis of hybrid improper ferroelectric s repre- 
sents a very promising route to the discovery of new 
magnetoelectric materials. The implications of rotation- 
driven ferroelectricity for the design of multifunctional 
materials were recently discussed in Ref. |75|. 

2.2.3. Inducing the linear ME ejfect in a topological in- 
sulator 

Although a^^ is much smaller than the spin counter- 
part for the conventional magnetoelectric s considered so 
far, Coh, et ah, argued that it may be possible to obtain 
very large orbital magnetoelectric couplings in materi- 
als close to a topological insulator (551 state. 6 is usu- 
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ally small in ordinary insulators, but it is equal to n in 
topological insulators (TI). This might lead one to think 
that TIs would exhibit a large linear ME, however the 
gapless surface states in TIs cancel the magnetoelectric 
contribution from the bulk and hence the total magneto- 
electric response is zero 158]. Note that it is time rever- 
sal symmetry that topologically protects these surface 
states, and equivalently, sets a to zero. However, one 
can imagine a material that does not have time reversal 
symmetry but is close to a TI state. Such a state can be 
obtained by breaking the time reversal symmetry of a 
TI by a small perturbation. This would lead to a large 
but nonquantized 6, namely - n. There are no topo- 
logically protected surface states in such a material so 
the linear ME response from the bulk can be observed. 
Coh and co-workers 1581 demonstrated the validity of 
this idea by breaking the time reversal symmetry in the 
well-studied TI Bi2Se3 by a fictitious Zeeman field and 
calculated from first principles. While this provides 
a proof of principle, discovering a material with strong 
orbital magnetoelectric response is currently one of the 
grand challenges in the field. 

3. Beyond Linear Magnetoelectric Effect: Spin - 
Phonon/Lattice Coupling 

The linear magnetoelectric eff'ect is the lowest or- 
der coupling between the electric and magnetic de- 
grees of freedom, however, it is not the only pos- 
sible coupling. Given the rather restrictive symme- 
try requirements of the linear eff'ect, researches have 
sought ways to achieve nontrivial magnetoelectric ef- 
fects from higher-order coupling of the electric and 
magnetic fields. One such higher-order coupling that 
has played a central role in the search for magneto- 
electrics in recent years is the biquadratic coupling, 
F ~ E^H^, which is allowed by symmetry in any mag- 
netic point group. 

One reason this term has been examined is because 
there is a simple physical origin of the term, which 
can be used as a guide to search for new materials and 
new eff'ects. The microscopic mechanism stems from 
the coupling of spins to the zone center, polar force- 
constants, often approximated in the literature by the 
following expression [i76il 

oj ^ copM-^y{Si • Sj), (11) 

where ojpm is the phonon frequency in the paramag- 
netic state well above a magnetic ordering temperature 
where local spins have formed but can be considered 
as uncorrelated, y is the spin-phonon coupling constant. 



and {Si ' Sj) is the nearest-neighbor spin-spin correla- 
tion function (note that in general the phonon frequency 
may also depend on correlation functions between fur- 
ther neighbors). This spin-phonon coupling exists in 
magnetic materials as a consequence of the fact that the 
magnetic exchange parameters depend on ionic posi- 
tions |76|. In this section, we will give a brief introduc- 
tion to the microscopic mechanism, discuss advances 
in our understanding of spin-phonon coupling in a few 
example materials and discuss the applicability of first- 
principles methods. Additionally, we will see how spin- 
lattice coupling, that is, the interdependence of mag- 
netic and structural orderings (which also stems from 
higher-order couplings) has been used to create novel 
magnetoelectric responses. 

3.1. Spin-Phonon and Spin Lattice coupling 

We begin with a few words regarding nomenclature. 
In principle, spin-phonon coupling refers to a dynamic 
eff'ect in which spin-correlations influence phonon fre- 
quencies and/or vice versa. Lawes and co-workers per- 
formed an elegant set of experiments demonstrating 
how such coupling leads to distinct dielectric constant 
vs temperature behavior in antiferromagnetic and ferro- 
magnetic materials f7T|. The manner in which mag- 
netic correlations influence the dielectric behavior of 
an insulator follows from examination of the expres- 
sion for the ionic contribution to the dielectric constant, 
Eq. [5j which goes like the inverse phonon frequency 
squared. Therefore, intrinsic magnetopermittivity has 
its origin in spin-phonon coupling. However, most stud- 
ies on spin-phonon coupling in the field of multiferroics 
are really discussing spin-lattice coupling: the eff'ect of 
the magnetic order/correlations on the force constant, 
an inherently static quantity. This misnomer carried 
over from the field of ferroelectrics, in particular the 
soft-mode theory of ferroelectricity in which a soft po- 
lar "phonon" is central to the theory (in reality, it is a 
softening of the lattice, that is, the force constant). We 
won't discuss these subtleties any further but unless oth- 
erwise noted we are discussing the coupling of spins to 
the force constants of the crystal. 

Techniques for the ab initio calculation of phonon 
frequencies and magnetic exchange parameters are well 
developed but researchers have only recently started cal- 
culating spin-phonon coupling ||78l |79l. The starting 
point has been the following model for the total energy 
of the system: 

^ = 2 ^ ^nmUnaUmp " ^ Jij^i ' + Eq, (12) 

nma/3 ij 
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where the first term is the internal elastic energy with 
Cnm as a bare force constant matrix, that is, the force 
constants in the paramagnetic state (note, this is not the 
same as the non-magnetic state, here we mean that lo- 
cal moments have formed but are uncorrelated) and Una 
denotes the displacement of atom n from an equilibrium 
position in the Cartesian direction a. The second term is 
the magnetic energy represented by a Heisenberg model 
(a good approximation for a large class of magnetic in- 
sulators) with Si as the spin of magnetic ion / and Jij as 
the exchange parameters. Finally, Eq is that part of the 
energy that doesn't depend on spins or atomic displace- 
ments. 

The idea (introduced decades ago |76|) which ac- 
counts for spin-phonon coupling is that in general the 
Ji/s not only depend on the positions of magnetic ions 
but also on the positions of the surrounding nonmag- 
netic atoms such, as the neighboring oxygen atoms. For 
small displacements, we can expand Jij in a Taylor se- 
ries up to second order in the displacements, for which 
it has been recently suggested 178] that the expansion 
should be done in the symmetrized basis functions, 77, 



'J 



(13) 



rather than in the ionic positions of two spins and an- 
gle w.r.t. the connecting anion. Here, n and m label ir- 
reducible representations and a and p label directions. 
J^. are exchange parameters for atoms in their equilib- 
rium positions. Here the assumption is that the first- 
order term is zero, which is true unless the crystal is 
a linear magnetoelectric (see Section [2]). By substitut- 
ing ( p3] ) into ^V2\ and integrating over the magnetic 
degrees of freedom one finds that the elastic energy 
is given by ^ 2 
force-constant matrix iTTSl 



nmap ^nm^nai^m/s with the rcnormalizcd 



(Si • S,-) 



(14) 



(15) 



where oj is the phonon frequency. Since the bare force 
constants are not accessible from a static theory such 
as density functional theory, Cnm cannot be directly cal- 
culated. Instead, Ref. |78| suggested that one should 
evaluate the force-constant matrix for several ordered, 
low-energy magnetic configurations (i.e., the left hand 
side of Eq. [14]). This allows one to obtain the deriva- 



tives 



dvnadvmp ^^^^^^ them to Eq. 
is known for the ordered system). 
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(since (S/ • S^) 
analogous to how 



the magnetic exchange parameters, Jij, are calculated 
using total-energy calculations. Spin-phonon coupling 
parameters can then be found as demonstrated for the 
spinel ZnCr204 fTSll . 

3.2. EuTiOi, and the idea of phase competition 

In most materials the magnetic energy scale is much 
smaller than the structural one, therefore the frequency 
change of a phonon due to spin-phonon coupling is ex- 
pected to be small. Recently, however, it was shown 
how this small spin-phonon coupling could not only 
produce a large/colossal magnetocapacitive eff'ect by 
mean of the structural softness idea discussed in section 
|2](the diff'erence being, f niguez discussed a bilinear cou- 
pling while here it is a higher order coupling) but also 
how it can be exploited to produce a large magnetoelec- 
tric response. It was realized that if a system is close to a 
ferroelectric phase transition, the spin-phonon coupling 
can produce radically diff'erent behavior by driving the 
system back and forth across the structural/magnetic 
boundary. Elaborating further, we know that in the soft- 
mode theory of ferroelectricity the polar ground state is 
associated with an unstable polar phonon. Figure [5ja). 
The idea is to take a material in which this polar phonon 
is additionally coupled to the spin system and tune the 
bare phonon frequency to near zero. Figure [5jb). In this 
case the structural ground state depends on whether the 
dominant magnetic interactions are antiferromagnetic 
Figure [Sjc), or ferromagnetic. Figure [Sj^d) (note the sign 
of the second derivative is system dependent). 

Hence, in the region where the bare phonon 
frequency is small, the competition between the 
antiferromagnetic-paraelectric and the ferromagnetic- 
ferroelectric states (or between the antiferromagnetic- 
ferroelectric and the ferromagnetic-paraelectric states 
if the sign of the second derivative is opposite) can 
lead to colossal magnetoelectric response. Addition- 
ally, the boundary of this phase competition region is 
bordered by a region of large/colossal magnetocapaci- 
tive responses due to the softening of the lattice. The 
encouraging aspect of this design approach is that the 
starting material is the typical magnetic insulator that is 
antiferromagnetic and paraelectric. 

An example of a material with spin-phonon coupling 
is EuTiOs, which is a PE insulator that crystallizes in 
the perovskite structure in bulk. The Eu^^ spins have G- 
type AFM order. Experiments on EuTiOs 1 80] showed 
that not only is there a suppression of the dielectric con- 
stant at the magnetic Neel temperature, but also that an 
external magnetic field of sufficient strength to align the 
magnetic moments causes a sizable increase in the di- 
electric constant. First-principles calculations (251 [HD 
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Figure 5: (a) Sketch of energy vs. normal mode displacement for 
stable (left) and unstable (right) phonon modes, (b) The phonon 
frequency squared can be separated into contributions from a bare 
phonon frequency (the frequency in the paraelectric state) and a mag- 
netic contribution that involves spin-spin correlations, (c) and (d) If 
the bare frequency is set to zero by an control parameter such as strain, 
the stability of the phonon mode depends on the sign of these correla- 
tions. 



verified that spin-phonon coupling, as suggested by the 
experimentaHsts, is indeed the reason for this behavior. 

In EuTiOs, the magnetic ion Eu is on the A site of 
the perovskite structure, unHke typical magnetic per- 
ovskites such as SrMnOs or SrCoOs . It is tempting to 
assume that the strength of exchange interactions does 
not depend on the ferroelectric mode, which mostly 
consists of B cation displacement. In order to explain 
the strong spin-phonon coupling that nevertheless ex- 
ists in this material, Akamatsu et al. 18211 , employed first 
principles calculations and showed that there exists a 
hybridization between the half-filled Eu / orbitals and 
the formally unoccupied Ti d orbitals. This observa- 
tion, combined with the fact that the bottom of the con- 
duction band consists dominantly of Ti d states, lead 
the authors to conclude that there exists an AFM su- 
perexchange interaction that is mediated via the Ti d or- 
bitals. The strength of this interaction depends strongly 
on the position of the Ti ion, and therefore the magnetic 
order is strongly coupled with the ferroelectric mode. 
This novel superexchange mechanism and the manner 
in which it couples with other lattice distortions, such 
as oxygen octahedron rotations, is a topic of active re- 
search, not only in the context of EuTiOs but also other 
rare earth - transition metal oxides |[83ll84ll. 
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Figure 6: (a) Sketch of soft mode frequency squared vs. biaxial strain 
for a material that becomes ferroelectric under strain, (b) An equiv- 
alent sketch for a magnetic material that has spin - phonon coupling 
and a lower soft mode frequency in the FM state. There is a ~ 0.2% 
wide strain range where the FM state has a ferroelectric instability but 
the AFM state does not. 



So EuTiOs has the desired characteristics: its ground 
state is antiferromagnertic and paraelectric yet parallel 
spins soften the polar phonon. But how can the system 
be "dialed" into the region of the phase diagram where 
the phase competition takes place? Epitaxial strain is 
widely used to stabilize new or meta- stable structural 
phases in transition metal oxides. In the case of per- 
ovskites it is well known that the polar soft mode is 
strongly sensitive to epitaxial strain and can even drive 
it from a paraelectric to a ferroelectric ground state. Fig- 
ure [6ja). A well known example is SrTiOs, which is a 
quantum paraelectric in bulk but has been experimen- 
tally demonstrated to display ferroelectricity at room 
temperature in thin film form under biaxial strain |85|. 
What has recently been appreciated is that in a mate- 
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rial with spin-lattice coupling, not only new structural 
phases but also new magnetic phases can be obtained 
by strain engineering. Figure [6jb) shows how the soft 
polar mode frequency varies with strain in an antiferro- 
magnetic material that has spin - phonon coupling and 
for which the system with parallel spins, i.e., FM, has 
a softer mode frequency than in the AFM state, such as 
EuTiOs. 

We point out that there is a wealth of information in 
Figure |6]) that is only accessible from first-principles. 
First, regardless of magnetic order, the soft mode is sta- 
ble for small values of strain, and unstable for large 
values. Next, as strain increases and the lattice soft- 
ens, there is a paraelectric region where the system with 
parallel spins has a near-zero polar phonon frequency 
while this phonon remains relatively hard in the ground 
state AFM. This is a region of large/colossal magneto- 
capacitance. Notice, however, the exact strain value at 
which the paraelectric to ferroelectric transition occurs 
depends on magnetic order. Therefore there is an inter- 
mediate strain range for which a ferroelectric instability 
exists only if the magnetic state is ferromagnetic. It is 
in this intermediate region where there is the possibility 
for phase competition between and antiferromagnetic- 
paraelectric and ferromagnetic-ferroelectric phases. If 
the ground state remains AFM-PE, it should be possible 
to apply a magnetic field, which would align the spins, 
thereby driving the system into the ferroelectric phase 
(or apply an electric field to induce an electrical polar- 
ization thereby driving the system into the FM phase), 
obtaining a large magnetoelectric response. The com- 
plete strain phase diagram of EuTiOs, calculated from 
first principles, is presented in Fig. [7] Note that there 
exists an intermediate region with large magnetoelectric 
response for both compressive and tensile strains, and 
the ground state is ferromagnetic - ferroelectric for large 
strain whether it is compressive or tensile. The reason 
is that biaxial strain favors a ferroelectric state, and the 
ferroelectric displacement of the Ti atom suppresses the 
AFM exchange interaction. This prediction has been re- 
cently verified for EuTiOs thin films grown on DyScOs 
substrate, which provides a 1.1% tensile strain |86|. 

The subsequent experimental realization of a 
ferroelectric-ferromagnetic state in epitaxially-strained 
thin films of EuTiOs was a breakthrough in multiferroic 
design and a proof-of-principle material 1861 (see also 
the article in this issue by Martin and Schlom). Several 
fundamental scientific challenges remain, however, be- 
fore materials useful for practical devices can be iden- 
tified. First, as discussed the present method relies on 
the ability to apply large epitaxial strains to a material. 
Such epitaxial strain engineering is readily possible to- 
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Figure 7: Properties of EuTiOs under biaxial strain. Ground state is 
AFM - PE for small values of strain, and FM - FE for large values. 
There exist intermediate regions with large ME effect for both tensile 
and compressive values of strain. 

day where the application of 1-2% biaxial tensile and 
compressive strains is routine. This method, however, 
requires growing a thin film coherently on a substrate of 
a larger/smaller lattice constant. The maximum thick- 
ness for which the film can be grown coherently with 
the substrate is determined as a balance between the 
volume strain energy of the film and the energetic cost 
of introducing misfit dislocations or cracks. For thick- 
nesses on the order of 50 nm for oxide films strained by 
1%, misfit dislocations are typically introduced, relax- 
ing the strain imposed across the interface. For a vari- 
ety of technological reasons, it would be beneficial to 
identify a control mechanism (a "knob" ) that can fea- 
sibly work for larger volumes of material. A second 
challenge concerns the temperature behavior of the ma- 
terial. Ideally, for practical applications in electronic, 
magnetic, or optical devices, the characteristic operat- 
ing temperature of the material should approach room 
temperature. In principle the novel mechanism of which 
EuTiOs is a realization could operate at room tempera- 
ture, however EuTiOs itself orders magnetically at 5.5 
K. It is therefore important to identify alternative mate- 
rials with higher magnetic ordering temperatures (above 
room temperature, preferably). As we have been dis- 
cussing in this review, first-principles methods are play- 
ing an important role in the search for such systems. 

33. Raising the temperature: f ^ d 

Recently, Lee and Rabe used first-principles calcula- 
tions to study the spin-phonon coupling in a series of 
perovskites SrMOs (M = V, Cr, Mn, Fe, Co) U3l In 
these materials, the magnetism is due to 3d electrons, 
which hybridize much more strongly with their envi- 
ronment than 4/ electrons. Exchange interactions in 
SrMOs are thus stronger than in the rare earth titanates 
(such as EuTiOs), leading to higher magnetic ordering 
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temperatures. In particular, for SrMnOs and SrCoOs, 
the transition temperatures are 233 K and 305 K, re- 
spectively. The calculations of Lee and Rabe indicated 
large spin-phonon coupling for the soft polar mode of 
SrMnOs and SrCoOs crystals |14|. Furthermore, for 
SrMnOs the polar mode became unstable on changing 
the magnetic configuration from AFM to FM, indicat- 
ing a FE transtion. In the case of SrCoOs, the system is 
metallic and the magnetic ground state is ferromagnetic. 
Spin-phonon coupling destabilizes a zone center mode 
(which would be polar if the system was not a metal) 
when G-type AFM ordering is imposed. 

Lee and Rabe also identified the microscopic mech- 
anisms responsible for the the spin-phonon coupling in 
SrMnOs and SrCoOs. For both materials, the mech- 
anism is directly related to the exchange interaction 
between spins. For SrMnOs, the ground state is G- 
type AFM and the Mn atoms interact through superex- 
change, which is mediated by the oxygen atoms. The 
superexchange mechanism leads to the AFM nearest- 
neighbor exchange parameter given by, 

(16) 



(a) 
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Jh is the Hund coupling, A is the energy difference be- 
tween the oxygen 2p and Mn 3J orbitals and t is the 
hopping integral between Mn nearest neighbors, which 
is proportional to the eff'ective wave function overlap be- 
tween the two Mn atoms through their shared oxygen. 
The hopping integral, and therefore the magnitude of 
Jnn, is maximized for the cubic PE structure where the 
Mn-O-Mn angle is 180° (see Fig. [Sja), left). In the FE 
structure, Mn atoms move toward one of the neighbor- 
ing oxygens (see Fig. [Sja), right) thereby decreasing 
the Mn-O-Mn angle, which leads to a decrease of 
In addition, the reduction of the bond length between 
the Mn atom and one of the oxygens enhances the hy- 
bridization between the Mn ?>d and O 2p orbitals, which 
increases A and decreases Jnn further. This reduction 
of the AFM exchange interaction in the FE state is re- 
sponsible for the spin-phonon coupling. Indeed, if the 
FM state is enforced then the energy cost for the FE dis- 
tortion decreases because it diminishes the energetically 
unfavorable exchange interaction. 

For metallic SrCoOs, the FM ground state ordering 
was shown to be a result of the double exchange inter- 
action |87|. Similarly, as in the case of AFM superex- 
change for SrMnOs, the Zener double exchange interac- 
tion, which depends sensitively on d-p hopping, is max- 
imized when the Co-O-Co angle is 180° and therefore 
it is diminished when the Co ion is displaced due to a 
phonon [see Fig. [SJb)]. This leads to the softening of 
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Figure 8: a) Schematic picture illustrating the origin of the spin- 
phonon coupling for the polar mode in SrMnOs. In the PE cubic 
state, the Mn-O-Mn angle is 180° and this maximizes the AFM su- 
perexchange interaction. The FE distortion reduces the the Mn-O-Mn 
angle and hence reduces the strength of the exchange interaction, b) 
As in a) but for SrCoOa. The overlap between Co- J and O-p orbitals 
determines the strength of Zener double exchange. The FE distortion 
decreases this overlap and hence the strength of the exchange inter- 
action, c) Illustration of the spin-phonon coupling in EuTiOs. The 
Ti-mediated superexchange between the Eu 4/ spins depends on the 
fourth power of the hopping amplitude tfd to the Ti d orbitals that are 
(partly) directed towards Eu ions. Displacement of the Ti ion from the 
center of the unit cell due to the polar phonon mode changes tfd and 
hence weakens the exchange interaction. 



this phonon mode as the G-type AFM ordering is en- 
forced d. 

Lee and Rabe also showed that both SrCoOs and 
SrMnOs undergo various phase transitions as a func- 
tion of strain ESUSgl. SrCoOs, a ferromagnetic metal, 
becomes an antiferromagnetic insulator under large 
enough tensile or compressive strain. SrMnOs, which is 
an antiferromagnetic paraelectric, becomes a ferromag- 
netic ferroelectric under large enough tensile or com- 
pressive strain. In addition to these phases, SrMnOs 
can be stabilized with various different antiferromag- 
netic states and has a magnetic ordering temperature of 
~ 100 K (See Fig. [9]). 

The transition of SrMnOs to a FE - FM state can be 
explained in the same way as EuTiOs : strain favors fer- 
roelectricity, which suppresses the AFM exchange in- 
teractions. Experimental observation of a similar spin - 
lattice coupling in perovskite Sri_;cBa;cMn03 |90 | sup- 
ports this explanation. Increasing the Ba concentra- 
tion in Sri_;cBa;cMn03 increases the lattice constant and 
drives the system ferroelectric for x > 0.45. The criti- 
cal temperature for ferroelectricity is larger than that of 
antiferromagnetism and the polarization decreases sig- 
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Figure 9: Properties of SrMnOs under biaxial strain. (Reproduced 
from Ref. | 88 1.) (a) Polarization of various phases, and the ground 
state (black), (b) Magnetic ordering temperature of the ground state, 
(c) c lattice constant. Copyright (2010) by the American Physical 
Society. 

nificantly upon magnetic ordering. The authors explain 
this reduction in polarization by the fact that the energy 
gain due to AFM order is largest when the Mn - O - Mn 
angle is 180°, so the magnetic ordering favors a nonpo- 
lar state. 

3.4. Spin-phonon coupling from first principles: which 
exchange -correlation functional to choose? 

First-principles calculations of magnetism in tran- 
sition metal oxides are usually performed using the 
T>¥T+U [91 1 method, where a Hubbard-type U term 
(partially) takes into account strong correlations be- 
tween localized electrons. Though practical, the ap- 
proach does have some serious drawbacks, including 
the fact that there is no obvious way to choose an ap- 
propriate value for U (there has been some progress 
on approaches in which U is evaluated self-consistently 
for the system of interest 1921 ). The standard pro- 
cedure usually involves varying U until some chosen 
quantity (such as the band gap, for example) matches 
the experimental value. Care must be taken to check 
the U dependence of the calculated properties of inter- 
est, as we will discuss further below. There are (be- 
yond DFT) techniques which explicitly capture the in- 



teraction between localized electrons, such as Dynam- 
ical Mean Field Theory (DMFT) |W, 94] or Quantum 
Monte Carlo (QMC) |i95J. DMFT in particular has 
provided many insights into the electronic structure of 
strongly correlated materials f96l|97l|98l|99l. However, 
as powerful as these techniques are, the calculation of 
forces is difficult, an obvious problem if the aim is to 
study spin-phonon coupling. 

Hong and co-workers | lOOJ studied the spin-phonon 
coupling in a series of perovskite transition metal oxides 
using both PBEsol-ht/ and the Heyd-Scuseria-Ernzerhof 
(HSE) hybrid functional 1 101 1. Hybrid functionals in- 
corporate a fraction of exact exchange from Hartree- 
Fock theory and are generally accepted to provide an 
improved description of the electronic structure of many 
materials; the Heyd-Scuseria-Ernzerhof (HSE) hybrid 
functional IIOII 11021 in particular has been shown to 
provide an improved description of the properties of 
transition metal oxides II103II (including phonon fre- 
quencies 11041 11051) with the advantage of being less 
computationally demanding than other hybrids, such as 
PBEO. However, the calculation of some quantities, 
such as phonon frequencies can still be prohibitively 
expensive. Hong, et ah, developed a procedure in 
which the value of U is determined by fitting to the 
HSE results for the diff'erence in total energy between 
diff'erent magnetic configurations. In testing their ap- 
proach, they found that the magnitude of spin-phonon 
coupling was very sensitive to the choice of U for cer- 
tain modes. They showed that the frequency shifts be- 
tween diff'erent magnetic arrangements calculated with 
PBEsol-ht/ (with the U fitted to the HSE data) were in 
very good agreement with the frequency shifts calcu- 
lated with HSE directly. This procedure therefore rep- 
resents a practical and convenient way to increase the 
reliability of such calculations without increasing the 
computational expense. 

4. The magnetoelectric effect in multiferroic het- 
erostructures and superlattices 

The discussion so far has been focused on homoge- 
neous materials and intrinsic mechanisms. Recent ad- 
vances in epitaxial growth techniques, however, have 
opened paths to create artificial multiferroic systems 
by forming layered heterostructures composed of mag- 
netic and dielectric materials. Such composite multifer- 
roics have attracted much attention, since they may ex- 
hibit a magnetoelectric (ME) eff'ect that is significantly 
larger than in single-phase multiferroic s. In general, 
we can distinguish between two types of ME coupling 
that occur in such heterostructures: (1) strain-mediated 
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ME effect, and (2) interfacial ME effect. In the for- 
mer type, the only effect of the interface is to strain the 
hetero structure constituents, which allows for control of 
their bulk properties. In this context, bilayers composed 
of piezoelectric and piezomagnetic materials have been 
intensively investigated. Here, the piezoelectric (piezo- 
magnetic) layer can be strained by application of an 
electric (magnetic) field. This strain is 'transmitted' to 
the neighboring piezomagnetic (piezoelectric) layer and 
changes its magnetization (polarization) leading to ME 
coupling I1Q6II1Q7||1Q8II . First-principles studies of this 
type of ME effect are usually concerned with the strain 
dependence of bulk properties of piezoelectric 11091 and 
piezomagnetic 1 1 10] materials and will not be discussed 
here. Instead, we focus on the second type of ME cou- 
pling, which is purely due to interfacial effects. 

Using first-principles calculations it was shown that 
a large ME effect is produced at the Fe/BaTiOs 
interface llllll I112L The Ti02 -terminated interface is 
stablel llT| and a magnetic moment opposite in direc- 
tion to the Fe magnetization is induced at the interfacial 
Ti atoms. The value of the Ti moment depends, how- 
ever, on the polarization of BaTiOs (BTO) being equal 
to 0.18 yUfi for the polarization pointing away from the 
interface and 0.40 //g for the polarization pointing to- 
ward the interface |111|. Hence, application of an elec- 
tric field that switches the polarization of the BTO leads 
to a change of interfacial magnetization. In order to 
quantify this effect, a surface ME coefficient was in- 
troduced, as = AMs/Ec, where AMg is the change of 
interface magnetization and Ec is the coercive field of 
the ferroelectric layer. Using the above results and the 
experimental value for the the coercive field of BTO, the 
surface ME coefficient was estimated to be 2 x 10"^^ G 
cm^/V [1111. This value, however, strongly depends on 
the oxygen stoichiometry at the interface 1 114|. 

The physical mechanism responsible for this inter- 
face ME effect can be understood by analyzing the 
orbital-resolved local density of states (DOS)| 1 1 1 1. Fig. 
[TOl shows that the Ti 3d band is centered at about 2 eV 
above the Fermi level and overlaps strongly with the mi- 
nority Fe 3d band, which has a significant weight in this 
energy region. This leads to hybridization between the 
Ti 3d and and minority Fe 3d states, which results in 
formation of bonding states that are pushed down in en- 
ergy and peaked just below the Fermi level (the peaks 
are indicated by arrows in Fig. [TOja) and[TOjb)). This 
results in a larger occupation of Ti minority states as 
compared to the majority states (the majority and mi- 
nority channels are determined by Fe), resulting in a 
magnetic moment on Ti that is antiparallel to the mag- 
netic moment of Fe BlllL If the polarization in BTO 
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Figure 10: Orbital-resolved DOS for interfacial atoms in a Fe/BaTiOs 
multilayer for net polarization of BaTiOa layer pointing toward (away 
from) the top (bottom) interface: (a) Ti 3d, (b) Fe 3d, and (c) O 2p. 
Majority- and minority-spin DOS are shown in the upper and lower 
panels, respectively. The solid and dashed curves correspond to the 
DOS of atoms at the top and bottom interfaces, respectively. The 
shaded plots are the DOS of atoms in the central monolayer of (b) 
Fe or (a),(c) Ti02 which can be regarded as bulk. The vertical line 
indicates the Fermi energy (Ef). From Ref. lllll . Copyright (2006) 
by the American Physical Society. 



points toward the interface the Fe-Ti bond length be- 
comes shorter, which enhances the hybridization and 
thus pushes the minority bonding state to even lower en- 
ergies so that it becomes more populated, resulting in a 
large magnetic moment on Ti. On the other hand, if the 
polarization points away from the interface, the Fe-Ti 
bond length increases, which weakens the hybridization 
and pushes the minority bonding state to higher ener- 
gies so that it is less populated and the Ti magnetic mo- 
ment decreases lllllll . First-principles calculations have 
revealed that such a 'hybridization-driven' ME effect 
also occurs for many other heterostructures, including 
Fe/PbTiOs (PTO) (III, Co2MnSi/BaTi03 |115| and 
Fe304/BTO ESI. 

Rondinelli et ah [111} proposed another mechanism 
to drive an interfacial ME effect based on the response 
of free electrons in a metal to an external electric field. 
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When a metal is placed into an electric field (E), the free 
electrons flow towards the surface and screen the elec- 
tric field over the screening length of the metal. In the 
case of a ferromagnetic metal this accumulated surface 
charge is spin polarized, which induces a net change in 
magnetization at the interface. Such an ME efl'ect driven 
by spin-dependent screening was demonstrated to occur 
at surfaces of ferromagnetic metals including Fe, Co, 
Nil. 118], SrRuOsLUaj), and halfmetallic (Fe304) 1 120 J . 
but the estimated ME coefl&cient is very small. The ME 
efl'ect may, however, be substantially enhanced at inter- 
faces between a metallic ferromagnet and an insulator 
due to the dielectric response of the insulator that effec- 
tively increases the electric field that must be screened 
by the ferromagnet. In particular, the induced surface 
charge, cr = sqsE, is proportional to the dielectric con- 
stant (s) of the insulator. Therefore, a large ME eff'ect 
can be expected at the interfaces between metallic ferro- 
magnets and high-6: dielectrics. Indeed, first principles 
calculations for the SrRuOs/SrTiOs (SRO/STO) inter- 
face in the presence of an electric field revealed a large 
ME response with ^ 2 x 10"^^ G cm^A^ ifTTTIl . The 
origin of this eff'ect was shown to be the spin-polarized 
screening which changes the occupation of the interfa- 
cial Ru-4J states, resulting in a change in local magnetic 
moment. |TTT| 

An even stronger ME eff'ect can be expected at the in- 
terface between ferromagnetic metals and ferroelectric s. 
In this case, changing the direction of polarization by 
application of an electric field changes the sign of po- 
larization charges at the interface and therefore modifies 
the induced surface charge in the ferromagnet at the in- 
terface. As the induced surface charge is spin-polarized, 
this induces a change in the interface magnetization. 
First-principles calculations show that this mechanism 
leads to the ME eff'ect at the SrRuOs/BaTiOs interface 
with as ^2x lO'^^ q ^^ijy 23 x IQ-^^ G cm^/V 
11191 , which is two orders of magnitude larger than for 
SrRuOs/SrTiOs. 

As discussed above, the screening of an electric field 
at metallic surfaces changes the carrier concentration in 
the surface region and thus can be thought of as electro- 
static doping. In particular, if we consider the surface 
of a metallic alloy then one can eff'ectively change the 
position of the system on the concentration phase dia- 
gram by applying an external field. This may lead to 
a dramatic eff'ect if the system sits close to the phase 
boundary, since then it would be possible to switch be- 
tween two phases using an electric field. In particular, 
a gigantic ME eff'ect can be achieved when the system 
sits at the boundary between FM and AFM phases. 

Lai_;,A;,Mn03 (LAMO) (A = Ca, Sr or Ba) alloys 



show a rich phase diagram as a function of doped car- 
rier concentration, which consists of diff'erent resistive, 
magnetic as well as structural phases. As an exam- 
ple we show a phase diagram for Lai_;cCa;cMn03 (see 
Fig. [TTJa)). For our purposes, the most interesting is 
the phase boundary between the FM metallic phase and 
the (AFM) insulating phase at concentrations of about 
0.5. Based on the above discussion, an external electric 
field may switch between these phases at the surface of 
Lao.sCao.sMnOs. Similarly, as in the case of induced 
surface magnetization, a larger eff'ect may be expected 
at the interface with an FE material since an electric 
field produced by polarization charges at the surface of 
a ferroelectric can be much larger than experimentally 
accessible electric fields. Indeed, first principles calcu- 
lations for (LAMO)/BTO interfaces have shown that 
switching the polarization of BTO by an electric field 
switches the magnetic order at the alloy side of the in- 
terface between AFM and FM (see Fig.[TT];b)) ifmll . 

Heterostructures built from metallic and ferroelectric 
materials provide a promising avenue of research, how- 
ever one must approach the first-principles treatment of 
such systems with extra care. The well-known 'band 
gap error' (the underestimation of the band gaps of 
semiconductors and insulators 11231 1124II ) may in par- 
ticular complicate calculations of interfaces. Stengel 
and co-workers have discussed in detail 1125 1 that the 
possible pathological alignment of the Fermi level of 
the metal with the conduction band of the ferroelectric 
at an interface could result in the incorrect prediction 
of a metallic heterostructure. The ferroelectric layer, 
which should be insulating, becomes metallic because 
its conduction band is occupied by free charges from 
the metal. Fundamental studies of the physics of such 
interfaces, as well as studies of the accuracy of diff'er- 
ent theoretical treatments (the accuracy of hybrid func- 
tional for band alignment, for example), are currently 
active areas of research. 

Recently, manipulation of magnetocrystalline 
anisotropy (MCA) by applying an electric field has 
been achieved for many FM metallic interfaces (sur- 
faces) |[T26l [127^. 128] as well as at the interface 
(surface) of dilute magnetic semiconductors 1129111301 . 
Interestingly, a significant change in perpendicular 
MCA energy (~ 40 %) by an applied voltage at the 
Fe/MgO interface 11271 and moreover the direct 
switching from out-of-plane to in-plane magnetiza- 
tion has been achieved in Au/Feo.8Coo.2/MgO |128|. 
Further, a modification of the out-of-plane magnetic 
anisotropy via piezoelectric response has been shown 
for the epitaxially grown magnetite films on the BaTiOa 
substrate lfT3ll . 
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Modulation of MCA by an electric field was also re- 
ported by ab initio studies | nH \lM [ESI [ml HSU 
11361 . In Fig. 12 ^a) we show a calculated electric field 
dependence of MCA for the Fe/MgO interface 11361 . 
In order to understand this eff'ect, we remind the reader 
again that MCA is a result of spin-orbit coupling. Us- 
ing second-order perturbation theory and neglecting the 
spin-flip terms, in the case of a fully occupied majority 
d spin channel (as in case of Fe), the MCA energy can 
be written as, 



-MCA 



= E[100]-E[001] = -^-(SALh-S,{L\\ (17) 



where ^ quantifies the strength of the coupling between 
L and S . Equation 
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' shows that the MCA energy varies 
with the orbital moment anisotropy (Am/). Indeed, a 
change of the MCA for the Fe/MgO interface is accom- 
panied by a change of the orbital moment anisotropy 



(see Fig.[T2[a)) 11361 . The origin of the change in the 
orbital moment anisotropy is found to be the relative 
change in the occupation of the d orbitals. As shown in 
Fig.p^b), when the applied electric field is away from 
the Fe layer, electrons want to move away from the in- 
terface region and therefore the occupation of the dxz 
and dyz orbitals, which have a large charge density at 
the interface, is reduced. At the same time, the occu- 
pation of the dxy orbitals (which have a smaller charge 
density at the interface) increases 11361 . This leads to 
a decrease in (L)^, which results in a reduction of Am/ 
and hence E^^^ (see Fig. 12 a)). A similar mechanism 
also explains the change in MCA at FM metal surfaces 
ifTTSl \T33i . On the other hand, for the Fe/BaTiOs in- 
terface it was shown that the change of the orbital mo- 
ment anisotropy (and therefore MCA) is due to modi- 
fication of the hybridization between Fe-3 J states and 
Ti-3 J states induced by the change of Fe-Ti bond length 
under polarization reversal 11321 . 



5. Outlook 

As the field moves forward, the focus of first- 
principles researchers working on the functional prop- 
erties of complex oxides is increasingly on the proper- 
ties of strongly correlated electron materials. We briefly 
touched on the technical difficulties associated with the 
first-principles treatment of such materials and men- 
tioned a couple of the advanced methods currently avail- 
able for calculations on strongly correlated electron sys- 
tems. The problem is that these advanced techniques are 
currently the province of a relatively small group of ex- 
perts and therefore need to become much more widely 
available and adopted if we are to bring the full force of 



the first-principles armory to bear on the study of mag- 
netoelectric phenomena in strongly correlated materials. 
Luckily, the resurgence of interest in the nickelates and 
cuprates has pushed these issues to the forefront of first- 
principles research on complex oxides. 
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6. Appendix 

6.1. Appendix I: The Dzyaloshinskii-Moriya (DM) in- 
teraction 

The focus of many studies on magnetoelectric ma- 
terials has been the the Dzyaloshinskii-Moriya (DM) 
interaction, also known as antisymmetric exchange (as 
opposed to symmetric Heisenberg exchange). The DM 
interaction between two spins S / and S j has the form 



£-- = 4..(^,x^,) 



(18) 



where is the DM vector. In 1957, Dzyaloshinskii 
pointed out 11371 that a free energy contribution of the 
form 

(19) 



F^^ = D.(MxL), 



where M is the magnetization, L is the antiferromag- 
netic vector and D is a material specific vector, is al- 
lowed by symmetry in certain materials. Shortly after 
this, Moriya 1 138] showed that the spin-orbit correction 
to Anderson superexchange leads to the microscopic in- 
teraction in Equation [TSj which in turn leads to the free 
energy contribution in Equation [T9j 

Since the DM interaction is a relativistic effect, it is 
usually weak in the sense that the energy scale it in- 
volves is much smaller than the energy scale of Heisen- 
berg exchange. However, it is still important because it 
gives rise to the weak ferromagnetism (wFM) observed 
in some antiferromagnets: any nonzero D causes the 
spins in an antiferromagnet to cant (typically by only 
a few degress) in order to minimize E^^ . This leads 
a small magnetic moment to appear in a material that 
would otherwise have zero macroscopic magnetization. 
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6.2. The linear magneto electric ejfect and weak- 
ferromagnetism 
Two interesting questions concerning how the elec- 
trical polarization couples to the magnetization have 
arisen in recent discussions of multiferroics. The first 
concerns the implications of a ferroelectric distortion 
inducing weak ferromagnetism 1631 while the second 
concerns ferroelectric distortion inducing linear mag- 
netoelectric eff'ect. What becomes apparent after a 
simple Landau exercise is that two diff'erent systems 
with identical magnetic point groups in a ferroelectric- 
magnetic phase can display qualitatively diff'erent cou- 
pling physics between the equilibrium magnetization 
and equilibrium polarization (as opposed to small 
changes in the magnetization and small changes to the 
polarization, i.e., the response to small fields, which 
must be the same since they have identical magnetic 
point groups). We will see how in one case an in- 
crease in the amplitude of the equilibrium polarization 
will lead to an increase in the magnitude of the equi- 
librium magnetization, but not in the linear magneto- 
electric eff'ect, while in a second case an increase in the 
equilibrium polarization leads to the exact opposite sit- 
uation. All of this becomes apparent by considering the 
relevant paraelectric reference structure and asking the 
question whether the PE structure allows a linear mag- 
netoelectric eff'ect or the eff'ect of weak- ferromagnetism, 
which in a paraelectric structure are mutually exclusive. 

6.2.1. Approach: the paraelectric reference structure 

In the present context, considering a paraelectric ref- 
erence structure is extremely illuminating. This is be- 
cause in a paraelectric phase, the linear magnetoelec- 
tric eff'ect and weak-ferromagnetism are mutually ex- 
clusive by symmetry I139l l3l l. To see this it is clearest 
to revert to Landau's original formulation of the linear 
ME eff'ect T{E, H) = -aijEiHj where E and H are the 
applied fields. It is obvious that if the magnetic point 
group contains either space inversion symmetry, J, or 
time-reversal symmetry, then the linear magnetoelec- 
tric eff'ect is zero. In a paraelectric phase, i.e., a phase 
containing a symmetry element that takes P ^ -P 
and thereby enforces P = 0, the only allowed sym- 
metry element that is compatible with the presence of 
a linear magnetoelectric eff'ect is the composite opera- 
tion of time-reversal R followed by space inversion, or 
RI. If the magnetic point group contains RI, it displays 
a linear magnetoelectric eff'ect, yet a net magnetization 
is forbidden as the system must remain invariant under 
RI\ RIM -M, therefore M =0. 

In a ferroelectric phase, where there is no required 
symmetry that enforces P -P, both the linear magne- 



toelectric eff'ect and weak-ferromagnetism are possible. 
A Landau argument, however, tells us nothing about the 
origin of these terms and hence nothing on the eff'ects 
of tuning. Yet, if one considers a paraelectric reference 
structure for the particular ferroelectric phase of inter- 
est, a lot can be learned. For example, if in the PE ref- 
erence structure the linear ME eff'ect is allowed, then 
the introduction of the ferroelectric distortion does not 
qualitatively change this, and therefore one doesn't ex- 
pect the linear ME eff'ect to scale linearly with the am- 
plitude of the inversion- symmetry breaking distortion. 
On the other hand, in this PE reference structure, weak- 
ferromagnetism is forbidden. The ferroelectric distor- 
tion then qualitatively changes the situation, and one 
does expect that as the size of the inversion-symmetry 
breaking distortion is increased, the size of the weak- 
ferromagnetism increases. The exact opposite conclu- 
sions are reached if the PE reference structure of inter- 
est displays weak-ferromagnetism, rather than the linear 
ME eff'ect. 

6.2.2. The result 

We wish to search for an antiferromagetic-paralectric 
reference structure having a group- subgroup relation to 
the FM-FE phase of interest. There are two cases that 
can arise. 

(1) The AFM-PE structure allows a linear magneto- 
electric term, ^, 



1 



1 



1 



TAiP, M) = -apP' + -a^M' + -bpR"" - ^L^PM (20) 

Since we are considering only the situation where a 
magnetization arises due to the coupling to some other 
order parameter, i.e., am > 0, we only need to consider 
0{M^). 

We now consider the value of the magnetization for 
a fixed value of the electrical polarization and find 
dTAldM = a^M-^LoP = 0^ M = ^LoP/a^, and then 
consider small changes of P due to an applied electric 
field, that is, F ^ Pq -\- AP 



M = Mo + AM 



with 



Mo = [— ) Po and AM = [— | AP (21) 



We see therefore that the equilibrium magnetization. Mo 
is proportional the equilibrium value of the electrical 
polarization 

Mo oc Po, 
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. The polarization vanishes in the paraelectric phase (if 
the paraelectric phase exists) and is directly proportional 
to the "size" of the ferroelectric distortion. Notice that if 
the direction of the equilibrium electrical polarization, 
Po, switches to its symmetry equivalent direction the 
equilibrium magnetization, Mq, switches 180° as well, 
and this information is not accessible from a Landau 
theory where the expansion is about the ground state 
ferroelectric phase. The second thing to notice is that a 
linear magnetoelectric effect, AM ~ AP, exists whether 
or not the polarization is nonzero. 

(2) In the second case of interest the AFM-PE struc- 
ture allows a weak-ferromagnetic term, J, 

rB{P,M) = ^apPH^amMH^bpP^^dLoM-^'LoP^M 

(22) 

We again consider the value of the magnetization for 
a fixed value of the electrical polarization and find 
OTbIOM = a^M+ JLo-f LqP^ = 0^ M = -JLo/a^- 
^^LoP^/am, and then considering the small changes of P 
due to an applied electric field (F ^ Pq -\- AP) we get 

^ M = Mo + AM 

with 

Mo = -^J^-^]pl and (23) 
AM = i^Po]AP 

We now see that the equilibrium magnetization, Mq is 
non-zero even in the case when the equilibrium value 
of the electrical polarization vanishes (that is, in the 
paraelectric phase, if the paraelectric phase exists). No- 
tice that if the direction of the electrical polarization 
switches to its symmetry equivalent direction the mag- 
netization does NOT switch. Again, this information is 
not accessible from a Landau theory where the expan- 
sion is about the ground state ferroelectric phase. The 
second thing to notice is that a linear magnetoelectric 
eff'ect becomes nonzero when the polarization switches 
on and increases in magnitude with Pq. 

References 

[1] M. L. Cohen, Predicting new solids and superconductors, Sci- 
ence 234 (1986) 549-553. 

[2] M. L. Cohen, Predicting useful materials. Science 261 (1993) 
307-308. 

[3] A. Franceschetti, A. Zunger, The inverse hand- structure prob- 
lem of finding an atomic configuration with given electronic 
properties, NATURE 402 (1999) 60-63. 



[4] N. A. Spaldin, W. E. Pickett, Computational design of mul- 
tifunctional materials. Journal of Solid State Chemistry 176 
(2003)615-632. 

[5] C. C. Fischer, K. J. Tibbetts, D. Morgan, G. Ceder, Predicting 
crystal structure by merging data mining with quantum me- 
chanics, Nat Mater 5 (2006) 641-646. 

[6] J. Hafner, C. Wolverton, G. Ceder, Toward computational ma- 
terials design: The impact of density functional theory on ma- 
terials research, MRS BULLETIN 31 (2006) 659-665. 

[7] R. Seshadri, S. L. Brock, A. Ramirez, M. Subramanian, M. E. 
Thompson, Advances in the development and growth of func- 
tional materials: Toward the paradigm of materials by design., 
MRS Bulletin 37 (2012) 682-690. 

[8] T. Kimura, T. Goto, H. Shintani, K. Ishizaka, T. Arima, 
Y. Tokura, Magnetic control of ferroelectric polarization, Na- 
ture 426 (2003) 55. 

[9] T. hisa Arima, Spin-driven ferroelectricity and magneto- 
electric eff'ects in frustrated magnetic systems. Journal of the 
Physical Society of Japan 80 (2011) 052001. 
[10] I. A. Sergienko, E. Dagotto, Role of the dzyaloshinskii-moriya 
interaction in multiferroic perovskites, Phys. Rev. B 73 (2006) 
094434. 

[11] S. Picozzi, K. Yamauchi, B. Sanyal, I. A. Sergienko, 
E. Dagotto, Dual nature of improper ferroelectricity in a mag- 
netoelectric multiferroic, Phys. Rev. Lett. 99 (2007) 227201. 

[12] I. A. Sergienko, C. Sen, E. Dagotto, Ferroelectricity in the 
magnetic E-phase of orthorhombic perovskites, Phys. Rev. 
Lett. 97 (2006) 227204. 

[13] H. Katsura, N. Nagaosa, A. V. Balatsky, Spin current and mag- 
netoelectric effect in noncoUinear magnets, Phys. Rev. Lett. 95 

(2005) 057205. 

[14] J. H. Lee, K. M. Rabe, Large spin-phonon coupling and 
magnetically induced phonon anisotropy in srmos perovskites 
(m = V,cr,mn,fe,co), Phys. Rev. B 84 (2011) 104440. 

[15] W. Zhong, D. Vanderbilt, K. M. Rabe, Phase transitions in 
batios from first principles, Phys. Rev. Lett. 73 (1994) 1861- 
1864. 

[16] W. Zhong, D. Vanderbilt, Eff'ect of quantum fluctuations on 
structural phase transitions in srtios and batios, Phys. Rev. B 
53 (1996) 5047-5050. 

[17] S. Prosandeev, L A. Kornev, L. Bellaiche, Magnetoelectricity 
in bifeos films: First-principles-based computations and phe- 
nomenology, Phys. Rev. B 83 (2011) 020102. 

[18] W. Ren, L. Bellaiche, Prediction of the magnetotoroidic ef- 
fect from atomistic simulations, Phys. Rev. Lett. 107 (2011) 
127202. 

[19] D. Sichuga, W. Ren, S. Prosandeev, L. Bellaiche, Chiral pat- 
terns of tilting of oxygen octahedra in zero-dimensional ferro- 
electrics and multiferroic s: A first principle-based study, Phys. 
Rev. Lett. 104 (2010) 207603. 

[20] W. F. Brown, R. M. Hornreich, S. Shtrikman, Upper bound 
on the magnetoelectric susceptibility, Phys. Rev. 168 (1968) 
574-577. 

[21] D. Vanderbilt, Private communication, ???? 

[22] R. E. Newnham, Phase Transformations in Smart Materials, 

Acta Crystallographica Section A 54 (1998) 729-737. 
[23] Y. Tokura, Critical features of colossal magnetoresistive man- 

ganites, Reports on Progress in Physics 69 (2006) 797. 
[24] Y. Tokura, Multiferroics - toward strong coupling between 

magnetization and polarization in a solid, JOURNAL OF 

MAGNETISM AND MAGNETIC MATERIALS 310 (2007) 

1145-1150. 

[25] C. J. Fennie, K. M. Rabe, Magnetic and electric phase control 
in epitaxial eutios from first principles, Phys. Rev. Lett. 97 

(2006) 267602. 



19 



[26] M. Fiebig, Revival of the magnetoelectric effect, Journal of 

Physics D: Applied Physics 38 (2005) R123. 
[27] M. Fiebig, N. A. Spaldin, Current trends of the magnetoelectric 

effect, EUROPEAN PHYSICAL JOURNAL B 71 (2009) 293- 

297. 

[28] R Curie, J. Physique 3 (1894) 393. 

[29] D. ASTROV, SOVIET PHYSICS JETP-USSR 11 (1960) 708- 
709. 

[30] L. Landau, L. Pitaevskii, E. Lifshitz, Electrodynamics of Con- 
tinuous Media, Butterworth-Heinemann, 1984. 

[31] A. S. Borovik-Romanov, H. Grimmer, Magnetic properties, 
volume D of International Tables for Crystallography, John 
Wiley & Sons, Ltd, pp. 105-149. 

[32] J. Iniguez, First-principles approach to lattice-mediated mag- 
netoelectric effects, Phys. Rev. Lett. 101 (2008) 117201. 

[33] W. COCHRAN, EFFECTIVE IONIC CHARGE IN CRYS- 
TALS, NATURE 191 (1961) 60-&. 

[34] P. Ghosez, J. Michenaud, X. Gonze, Dynamical atomic 
charges: The case of ABO(3) compounds, PHYSICAL RE- 
VIEW B 58 (1998) 6224-6240. 

[35] J. AXE, APPARENT IONIC CHARGES AND VIBRA- 
TIONAL EIGENMODES OF BATI03 AND OTHER PER- 
OVSKITES, PHYSICAL REVIEW 157 (1967) 429-&. 

[36] W. COCHRAN, R. COWLEY, DIELECTRIC CONSTANTS 
AND LATTICE VIBRATIONS, JOURNAL OF PHYSICS 
AND CHEMISTRY OF SOLIDS 23 (1962) 447-&. 

[37] A. Scaramucci, E. Bousquet, M. Mostovoy, N. A. Spaldin, 
Linear Magnetoelectric Effect by Orbital Magnetism, ArXiv 
e-prints (2012). 

[38] J. C. Wojdel, J. ffiiguez, Magnetoelectric response of multifer- 
roic bifeo3 and related materials from first-principles calcula- 
tions, Phys. Rev. Lett. 103 (2009) 267205. 

[39] J. C. Wojdel, J. Iiiiguez, Ab Initio indications for giant magne- 
toelectric effects driven by structural softness, Phys. Rev. Lett. 
105 (2010) 037208. 

[40] E. Bousquet, N. A. Spaldin, K. T. Delaney, Unexpectedly large 
electronic contribution to linear magnetoelectricity, Phys. Rev. 
Lett. 106 (2011) 107202. 

[41] G. Kresse, J. Furthmiiller, Efficient iterative schemes for ab 
initio total-energy calculations using a plane-wave basis set, 
Phys. Rev. B 54 (1996) 11169-11186. 

[42] K. Yosida, Theory of Magnetism, Springer, 2001. 

[43] Claude cohen tannoudji, B. Diu, F. Laloe, Quantum Mechanics, 
Wiley-Interscience, 1992. 

[44] M. Mostovoy, A. Scaramucci, N. A. Spaldin, K. T. Delaney, 
Temperature-dependent magnetoelectric effect from first prin- 
ciples, Phys. Rev. Lett. 105 (2010) 087202. 

[45] X.-L. Qi, T. L. Hughes, S.-C. Zhang, Topological field theory 
of time-reversal invariant insulators, Phys. Rev. B 78 (2008) 
195424. 

[46] A. M. Essin, J. E. Moore, D. Vanderbilt, Magnetoelectric po- 

larizability and axion electrodynamics in crystalline insulators, 

Phys. Rev. Lett. 102 (2009) 146805. 
[47] A. M. Essin, A. M. Turner, J. E. Moore, D. Vanderbilt, Orbital 

magnetoelectric coupling in band insulators, Phys. Rev. B 81 

(2010) 205104. 

[48] A. Malashevich, I. Souza, S. Coh, D. Vanderbilt, Theory of 
orbital magnetoelectric response. New Journal of Physics 12 

(2010) 053032. 

[49] T. Thonhauser, International Journal of Modem Physics B 25 

(2011) 1429-1458. 

[50] T. Thonhauser, D. Ceresoh, D. Vanderbilt, R. Resta, Orbital 
magnetization in periodic insulators, Phys. Rev. Lett. 95 (2005) 
137205. 

[51] D. Ceresoli, T. Thonhauser, D. Vanderbilt, R. Resta, Or- 



bital magnetization in crystalline solids: Multi-band insulators, 
chem insulators, and metals, Phys. Rev. B 74 (2006) 024408. 
[52] R. Resta, D. Ceresoli, T. Thonhauser, D. Vanderbilt, Orbital 
magnetization in extended systems, ChemPhysChem 6 (2005) 
1815-1819. 

[53] F. Wilczek, Two applications of axion electrodynamics, Phys. 

Rev. Lett. 58 (1987) 1799-1802. 
[54] M. Franz, High-energy physics in a new guise. Physics 1 

(2008) 36. 

[55] M. Z. Hasan, C. L. Kane, Colloquium : Topological insulators. 

Rev. Mod. Phys. 82 (2010) 3045-3067. 
[56] G. H. Wannier, The structure of electronic excitation levels in 

insulating crystals, Phys. Rev. 52 (1937) 191-197. 
[57] N. Marzari, D. Vanderbilt, Maximally localized generalized 

wannier functions for composite energy bands, Phys. Rev. B 

56 (1997) 12847-12865. 
[58] S. Coh, D. Vanderbilt, A. Malashevich, I. Souza, Chern-simons 

orbital magnetoelectric coupling in generic insulators, Phys. 

Rev. B 83 (2011)085108. 
[59] A. Malashevich, S. Coh, I. Souza, D. Vanderbilt, To be pub- 

hshed, 2012. 

[60] S. Bhattacharjee, E. Bousquet, P. Ghosez, Engineering multi- 
ferroism in canmos, Phys. Rev. Lett. 102 (2009) 117602. 

[61] E. Bousquet, N. Spaldin, Induced magnetoelectric response in 
puma perovskites, Phys. Rev. Lett. 107 (2011) 197603. 

[62] A. Filippetti, N. A. Hill, Coexistence of magnetism and ferro- 
electricity in perovskites, Phys. Rev. B 65 (2002) 195120. 

[63] C. J. Fennie, Ferroelectrically induced weak ferromagnetism 
by design, Phys. Rev. Lett. 100 (2008) 167203. 

[64] C. Ederer, C. J. Fennie, Electric- filed switchable magnetiza- 
tion via the dzyaloshinskii-moriya interaction: FeTiOs versus 
BiFeOs, J. Phys.: Condens. Matter 20 (2008) 434219. 

[65] D. L. Fox, J. R Scott, J. Phys. C: Solid State Phys. 10 (1977) 
L329-331. 

[66] T. Varga, A. Kumar, E. Vlahos, S. Denev, M. Park, S. Hong, 
T. Sanehira, Y. Wang, C. J. Fennie, S. K. Streiffer, X. Ke, 
P. Schiffer, V. Gopalan, J. F. Mitchell, Coexistence of weak fer- 
romagnetism and ferroelectricity in the high pressure LiNbOs- 
type phase of fetios, Phys. Rev. Lett. 103 (2009) 047601. 

[67] Y. Inaguma, M. Yoshida, T. Katsumata, J. Am. Chem. Soc. 130 
(2008) 6704-6705. 

[68] J. Y. Son, G. Lee, M. H. Jo, H. Kim, H. M. Jung, Y. H. Shin, J. 
Am. Chem. Soc. 131 (2009) 8386-8387. 

[69] A. Aimi, T. Katsumata, D. Mori, D. Fu, M. Itoh, T. Kyomen, 
K. Hiraki, T. Takahashi, Y. Inaguma, Inorg. Chem. 50 (2011) 
6392-6398. 

[70] N. A. Benedek, C. J. Fennie, Hybrid improper ferroelectric- 
ity: a mechanism for controllable polarization-magnetization 
coupling, Phys. Rev. Lett. 106 (2011) 107204. 

[71] A. B. Harris, Symmetry analysis for the Ruddlesden-Popper 
systems Ca3Mn207 and Ca3Ti207, PHYSICAL REVIEW B 
84 (2011). 

[72] E. Bousquet, M. Dawber, N. Stucki, C. Lichtensteiger, P. Her- 
met, S. Gariglio, J. Triscone, P. Ghosez, Improper ferroelec- 
tricity in perovskite oxide artificial superlattices. Nature 452 
(2008) 732-737. 

[73] J. M. Rondinelli, C. J. Fennie, Octahedral rotation-induced 
ferroelectricity in cation ordered perovskites. Adv. Mater. 24 
(2012) 1961. 

[74] M. W. Lufaso, P. M. Woodward, Prediction of the crystal struc- 
tures of perovskites using the software program SPuDS, Acta 
Cryst. B 57 (2001) 725-738. 

[75] N. A. Benedek, A. T. Mulder, C. J. Fennie, Po- 
lar octahedral rotations: a path to new multifunctional 
materials, J. Sohd State Chemistry In Press (2012) 



20 



http://dx.doi.Org/10.1016/j.jssc.2012.04.012. 

[76] W. Baltensperger, Influence of magnetic order on conduction 
electrons and phonons in magnetic semiconductors, Journal of 
Applied Physics 41 (1970) 1052-1054. 

[77] G. Lawes, A. Ramirez, C. Varma, M. Subramanian, Magne- 
todielectric eff'ects from spin fluctuations in isostructural fer- 
romagnetic and antiferromagnetic systems, PHYSICAL RE- 
VIEW LETTERS 91 (2003). 

[78] C. J. Fennie, K. M. Rabe, Magnetically Induced Phonon 
Anisotropy in ZnCr204 from First Principles, Phys. Rev. Lett. 
96 (2006) 205505. 

[79] R. Sabiryanov, S. Jaswal, Magnons and magnon-phonon in- 
teractions in iron, PHYSICAL REVIEW LETTERS 83 (1999) 
2062-2064. 

[80] T. Katsufuji, H. Takagi, Coupling between magnetism and di- 
electric properties in quantum paraelectric eutios, Phys. Rev. 
B 64 (2001) 054415. 

[81] R. Ranjan, H. S. Nabi, R. Pentcheva, Electronic structure and 
magnetism of eutio 3 : a first-principles study, Journal of 
Physics: Condensed Matter 19 (2007) 406217. 

[82] H. Akamatsu, Y. Kumagai, F. Oba, K. Fujita, H. Murakami, 
K. Tanaka, I. Tanaka, Antiferromagnetic superexchange via 
3d states of titanium in eutios as seen from hybrid hartree- 
fock density functional calculations, Phys. Rev. B 83 (2011) 
214421. 

[83] H. Akamatsu, K. Fujita, H. Hayashi, T. Kawamoto, Y. Kuma- 
gai, Y. Zong, K. Iwata, F. Oba, I. Tanaka, K. Tanaka, Crystal 
and Electronic Structure and Magnetic Properties of Divalent 
Europium Perovskite Oxides EuMOs (M = Ti, Zr, and Hf): 
Experimental and First-Principles Approaches, Inorg. Chem. 
51 (2012) 4560-4567. 

[84] T. Birol, C. J. Fennie, In preparation (2012). 

[85] J. H. Haeni, P. Irvin, W. Chang, R. Uecker, R Reiche, Y. L. Li, 
S. Choudhury, W. Tian, M. E. Hawley, B. Craigo, A. K. Tagant- 
sev, X. Q. Pan, S. K. Streiffer, L. Q. Chen, S. W. Kirchoefer, 
J. Levy, D. G. Schlom, Room- temperature ferroelectricity in 
strained srtio3, Nature 430 (2004) 758-761. 

[86] J. H. Lee, L. Fang, E. Vlahos, X. Ke, Y. W. Jung, L. E 
Kourkoutis, J.-W. Kim, P. J. Ryan, T. Heeg, M. Roeckerath, 
V. Goian, M. Bernhagen, R. Uecker, P. C. Hammel, K. M. 
Rabe, S. Kamba, J. Schubert, J. W. Freeland, D. A. MuUer, 

C. J. Fennie, P. Schiff'er, V. Gopalan, E. Johnston-Halperin, 

D. G. Schlom, A strong ferroelectric ferromagnet created by 
means of spin-lattice couphng. Nature 466 (2010) 954-958. 

[87] C. Zener, Interaction between the d shells in the transition 

metals, Phys. Rev. 81 (1951) 440^44. 
[88] J. H. Lee, K. M. Rabe, Epitaxial- strain-induced multiferroicity 

in srmno3 from first principles, Phys. Rev. Lett. 104 (2010) 

207204. 

[89] J. H. Lee, K. M. Rabe, Coupled magnetic -ferroelectric metal- 
insulator transition in epitaxially strained srcoos from first 
principles, Phys. Rev. Lett. 107 (2011) 067601. 

[90] H. Sakai, J. Fujioka, T. Fukuda, D. Okuyama, D. Hashizume, 
F. Kagawa, H. Nakao, Y. Murakami, T. Arima, A. Q. R. Baron, 
Y. Taguchi, Y. Tokura, Phys. Rev. Lett. 107 (2011) 137601. 

[91] V. I. Anisimov, F. Aryasetiawan, A. I. Lichtenstein, First- 
principles calculations of the electronic structure and spectra 
of strongly correlated systems: the Ida -i- u method. Journal of 
Physics: Condensed Matter 9 (1997) 767. 

[92] H. J. Kulik, M. Cococcioni, D. A. Scherhs, N. Marzari, Den- 
sity functional theory intranstion-metal chemistry: A self- 
consistent hubbard-u approach, Phys. Rev. Lett. 97 (2006) 
103001. 

[93] A. Georges, G. Kotliar, W. Krauth, M. J. Rozenberg, Dynam- 
ical mean-field theory of strongly correlated fermion systems 



and the limit of infinite dimensions. Rev. Mod. Phys. 68 (1996) 
12-125. 

[94] G. Kothar, S. Y. Savrasov, K. Haule, V. S. Oudovenko, O. Par- 
collet, C. A. Marianetti, Electronic structure calculations with 
dynamical mean-field theory. Rev. Mod. Phys. 78 (2006) 865- 
951. 

[95] W. M. C. Foulkes, L. Mitas, R. J. Needs, G. Rajagopal, Quan- 
tum monte carlo simulations of solids. Rev. Mod. Phys. 73 
(2001) 33-83. 

[96] E. Pavarini, S. Biermann, A. Poteryaev, A. I. Lichtensteni, 
A. Georges, O. K. Andersen, Mott transition and supression of 
orbtial fluctuations in orthorhombic 3d(l) perovskites, Phys. 
Rev. Lett. 92 (2004) 176403. 

[97] T. V. Ramakrishnan, H. R. Krishnamurthy, S. R. Hassan, G. V. 
Pai, Theory of insulator metal transition and cololssal magne- 
toresistance in doped manganites, Phys. Rev. Lett. 92 (2004) 
157203. 

[98] S. Biermann, A. Poteryaev, A. I. Lichtenstein, A. Georges, Dy- 
namical singlets and correlation-assisted peierls transition in 
VO2, Phys. Rev. Lett. 94 (2005) 026404. 

[99] M. J. Han, X. Wang, C. A. Marianetti, A. J. Milhs, Dynamical 
mean-field theory of nickelate superlattices, Phys. Rev. Lett. 
107(2011)206804. 
[100] J. Hong, A. Stroppa, J. Iniguez, S. Picozzi, D. Vanderbilt, 
Spin-phonon coupling efl'ects in transition- metal perovskites: 
A dft M and hybrid-functional study, Phys. Rev. B 85 (2012) 
054417. 

[101] J. Heyd, G. E. Scuseria, M. Emzerhof, Erratum: "hybrid func- 
tionals based on a screened coulomb potential" [j. chem. phys. 
[bold 118], 8207 (2003)], The Journal of Chemical Physics 
124 (2006) 219906. 

[102] T. M. Henderson, J. Paier, G. E. Scuseria, Accurate treatment 
of solids with the hse screened hybrid, Phys. Status Sohdi B 
248 (2011)767-774. 

[103] T. Archer, C. D. Pemmaraju, S. Sanvito, C. Franchini, J. He, 
A. Filippetti, P. Delugas, D. Puggioni, V. Fiorentini, R. Ti- 
wari, P. Majumdar, Exchange interactions and magnetic phases 
of transition metal oxides: Benchmarking advanced ab initio 
methods, Phys. Rev. B 84 (2011) 115114. 

[104] R. Wahl, D. Vogtenhuber, G. Kresse, SrTiOs and BaTiOs re- 
visited using the projector augmented wave method: Perfor- 
mance of hybrid and semilocal functionals, Phys. Rev. B 78 
(2008) 104116. 

[105] C. Franchini, A. Sanna, M. Marsman, G. Kresse, Structural, 
vibrational and quasiparticle properties of the Peierls semi- 
conductor BaBiOs: A hybrid functional and self-consistent 
GW-Fvertex-corrections study, Phys. Rev. B 81 (2010) 085213. 

[106] J. van Suchtelen, Product properties: a new application of com- 
posite materials. Philips Research Reports 27 (1972) 28. 

[107] L. Mitoseriuab, V. Buscagliac, Intrinsic/extrinsic interplay 
contributions to the functional properties of ferroelectric- 
magnetic composites, Phase Transitions 79 (2006) 1095. 

[108] C.-G. Stefanita, From Bulk to Nano, volume 117, Springer, 
2008. 

[109] R-E. Janolin, N. A. Pertsev, D. Sichuga, L. Bellaiche, Giant di- 
rect magnetoelectric effect in strained multiferroic heterostruc- 
tures, Phys. Rev. B 85 (2012) 140401. 

[110] P. Lukashev, R. F. Sabirianov, K. Belashchenko, Theory of the 
piezomagnetic effect in mn-based antiperovskites, Phys. Rev. 
B 78 (2008) 184414. 

[Ill] C.-G. Duan, S. S. Jaswal, E. Y. Tsymbal, Predicted magneto- 
electric efl"ect in Fe/batios multilayers: Ferroelectric control of 
magnetism, Phys. Rev. Lett. 97 (2006) 047201. 

[112] M. Fechner, I. V. Maznichenko, S. Ostanin, A. Ernst, J. Henk, 
P. Bruno, I. Mertig, Magnetic phase transition in two-phase 



21 



multiferroics predicted from first principles, Phys. Rev. B 78 
(2008) 212406. 

113] M. Fechner, S. Ostanin, I. Mertig, Effect of the surface polar- 
ization in polar perovskites studied from first principles, Phys. 
Rev. B 77 (2008) 094112. 

1 14] M. Fechner, S. Ostanin, I. Mertig, Eff'ect of oxidation of the ul- 
trathin fe electrode material on the strength of magnetoelectric 
coupling in composite multiferroics, Phys. Rev. B 80 (2009) 
094405. 

115] K. Yamauchi, B. Sanyal, S. Picozzi, Interface eff'ects at a 
half-metal/ferroelectric junction, Appl. Phys. Lett. 91 (2007) 
062506. 

116] M. K. Niranjan, J. P. Velev, C.-G. Duan, S. S. Jaswal, E. Y. 
Tsymbal, Magnetoelectric effect at the fe304/batio3 (001) 
interface: A first-principles study, Phys. Rev. B 78 (2008) 
104405. 

117] J. M. RondineUi, M. Stengel, N. A. Spaldin, Carrier-mediated 
magnetoelectricity in complex oxide heterostructures. Nature 
Nanotechnology 3 (2008) 46. 

118] C.-G. Duan, J. R Velev, R. F. Sabirianov, Z. Zhu, J. Chu, S. S. 
Jaswal, E. Y. Tsymbal, Surface magnetoelectric effect in ferro- 
magnetic metal films, Phys. Rev. Lett. 101 (2008) 137201. 

119] M. K. Niranjan, J. D. Burton, J. R Velev, S. S. Jaswal, E. Y. 
Tsymbal, Magnetoelectric effect at the srruos/batios (001) in- 
terface: An ab initio study, Appl. Phys. Lett. 95 (2009) 052501. 

120] C.-G. Duan, C.-W. Nan, S. S. Jaswal, E. Y. Tsymbal, Uni- 
versality of the surface magnetoelectric effect in half-metals, 
Phys. Rev. B 79 (2009) 140403. 

121] J. D. Burton, E. Y. Tsymbal, Prediction of electrically in- 
duced magnetic reconstruction at the manganite/ferroelectric 
interface, Phys. Rev. B 80 (2009) 174406. 

122] T. V. Ramakrishnan, Modelling colossal magnetoresistance 
manganites, J. Phys.: Condens. Matter 19 (2007) 125211. 

123] J. P. Perdew, M. Levy, Physical content of the exact kohn- 
sham orbital energies: Band gaps and derivative discontinu- 
ities, Phys. Rev. Lett. 51 (1983) 1884-1887. 

124] L. J. Sham, M. Schliiter, Density-functional theory of the en- 
ergy gap, Phys. Rev. Lett. 51 (1983) 1888-1891. 

125] M. Stengel, P. Aguado-Puente, N. A. Spaldin, J. Junquera, 
Band alignment at metal/ferroelectric interfaces: Insights and 
artifacts from first principles, Phys. Rev. B 83 (2011) 235112. 

126] M. Weisheit, S. Fahler, A. Marty, Y. Souche, C. Poinsignon, 
D. Givord, Electric field-induced modification of magnetism 
in thin-film ferromagnets. Science 315 (2007) 349. 

127] T. Maruyama, Y. Shiota, T. Nozaki, K. Ohta, N. Toda, 
M. Mizuguchi, A. A. Tulapurkar, T. Shinjo, M. Shiraishi, 
S. Mizukami, Y. Ando, Y. Suzuki, Large voltage-induced mag- 
netic anisotropy change in a few atomic layers of iron. Nature 
Nanotechnology 4 (2009) 158. 

128] Y. Shiota, T. Maruyama, T. Nozaki, M. Shinjyo, T. andShi- 
raishi, Y. Suzuki, Voltage-assisted magnetization switching in 
ultrathin fe80CO20 alloy layers, Appl. Phys. Express 2 (2009) 
063001. 

129] D. Chiba, M. Sawicki, Y. Nishitani, Y. Nakatani, F. Matsukura, 
H. Ohno, Magnetization vector manipulation by electric fields. 
Nature 455 (2008)515. 

130] K. Ohta, T. Maruyama, T. Nozaki, M. Shiraishi, T. Shinjo, 
Y. Suzuki, S.-S. Ha, C.-Y. You, W. Van Roy, Voltage con- 
trol of in-plane magnetic anisotropy in ultrathin fe/n-gaas(001) 
schottky junctions, Appl. Phys. Lett. 94 (2009) 032501. 

131] C. A. F Vaz, J. Hoffman, A.-B. Posadas, C. H. Ahn, Mag- 
netic anisotropy modulation of magnetite in fe3O4/batio3(100) 
epitaxial structures, Appl. Phys. Lett. 94 (2009) 022504. 

132] C.-G. Duan, J. R Valev, R. F Sabirianov, W. N. Mei, S. S. 
Jaswal, E. Tsymbal, Tailoring magnetic anisotropy at the fer- 



romagnetic/ferroelectric interface, Appl. Phys. Lett. 92 (2008) 
122905. 

[133] K. Nakamura, R. Shimabukuro, Y. Fujiwara, T. Akiyama, 
T. Ito, A. J. Freeman, Giant modification of the magnetocrys- 
talline anisotropy in transition-metal monolayers by an exter- 
nal electric field, Phys. Rev. Lett. 102 (2009) 187201. 

[134] M. Tsujikawa, T. Oda, Finite electric field effects in the large 
perpendicular magnetic anisotropy surface Pt/Fe/Pt(001): A 
first-principles study, Phys. Rev. Lett. 102 (2009) 247203. 

[135] H. Zhang, M. Ritcher, K. Koepemik, I. Opahle, F. Tas- 
nadi, H. Eschrig, Electric-field control of surface magnetic 
anisotropy: a density functional approach. New J. Phys. 11 
(2009) 043007. 

[136] M. K. Niranjan, C.-G. Duan, S. S. Jaswal, E. Y. Tsymbal, Elec- 
tric field effect on magnetization at the fe/mgo(001) interface, 
Appl. Phys. Lett. 96 (2010) 222504. 

[137] I. E. Dzyaloshinskii, Thermodynamic theory of weak ferro- 
magnetism in antiferromagnetic substances, Sov. Phys. - JETP 
5 (1957) 1259-1272. 

[138] T. Moriya, Anisotropic superexchange interaction and weak 
ferromagnetism, Phys. Rev. 120 (1960) 91. 

[139] E. A. Turov, Can the magnetoelectric effect coexist with weak 
piezomagnetism and ferromagnetism?, Physics-Uspekhi 37 
(1994) 303. 



22 




0.0 0.2 0.4 0.6 0.8 1.0 



Cax 




A. 




E (V/nm) E 



Figure 12: (a) Magnetic anisotropy energy (MCA, squares) and orbital 
moment anisotropy (Am/, triangles) of the Fe/MgO(001) interface as 
a function of electric field in MgO. (b) Induced charge density Ap = 
p{E) -p(0), in units of e /A? , at the interfacial Fe atom for E= 1 .0 V/nm 
in the x-z (010) plane (top panel) and the x-y (001) plane (bottom 
panel). Reprinted with permission from Ref. |136 |. Copyright (2010) 
American Institute of Physics. 



Figure 11: (a) Phase diagram of Lai-j^CajtMnOs as a function of dop- 
ing concentration x and temperature T. From Ref.| 122|. (b) Electri- 
cally induced magnetic reconstruction at the LAMO/BTO interface; 
the predicted change in order of the Mn magnetic moments (small ar- 
rows) from FM to A-type antiferromagnetic (A2) as the ferroelectric 
polarization (large arrows) of the BTO is reversed. From Ref. L 1211 . 
Copyright (2009) by the American Physical Society. 
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